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1.0 SUMMARY 


This report describes the results of Phase II of the V/STOL Pro- 
pulsion Control Analysis Program for the development of propulsion 
control technologies for achieving integrated V/STOL aircraft-engine 
controls. The Phase I study was focused on the development of a long- 
range technology plan, the selection of a representative baseline 
V/STOL propulsion system, and the development of propulsion control 
design concepts for the vertical flight regime. The current Phase 
II program retained the Phase I baseline propulsion system and 
developed propulsion control design concepts for the flight transition 
regime between vertical and horizontal flight. 

The General Electric HALS (Remote Augmented Lift System) concept 
was retained as a baseline engine and was used to establish typical 
operating conditions along Vertlcal-to-Horizontal (VTH) and Horizontal- 
to-Vertical (HTV) transition processes. Analytical studies were con- 
ducted at each selected operating condition to establish steady-state 
operating requirements, to define a recommended closed-loop regulator 
configuration, and to establish regulator gains for achieving stable 
control operation. Regression techniques were utilized to develop one 
set of gain schedules for the VTH transition and, a second set of 
schedules for the HTV transition. 

Baseline engine component performance characteristics were utilized 
to develop a component -level regression model of the baseline engine. 
Internal model aerothermo relationships were based on approximated real 
variable gas properties. The resulting non-linear transient engine 
model was combined with non-linear transient models of the actuators, 
augmentors, and nozzles to provide a non-linear simulation capability 
for evaluating the regulator gain schedules. 

Additional mathematical models were developed for the inlet and 
nozzles and for a typical V/STOL aircraft configured for RALS. The 
inlet and nozzle effects model represents inlet distortion effects on 
the fan and compressor and nozzle deflection effects on flow coefficient 


The aircraft model was combined with a simple pilot model and used to 
examine typical transition trajectories and their corresponding pro- 
pulsion control requirements* 

The long-range propulsion control technology plah was refined to 
focus on those technology requirements essential for achieving a piloted 
simulation capability at NASA-Ames in the 1984-85 time period. 

2.0 INTRODUCTION 


A V/STOL propulsion system must provide vertical thrust for takeoff 
and landing, horizontal thrust for conventional flight, and thrust com- 
ponents in each direction for achieving a controlled transition between 
the two flight regimes. Concurrent flight path and attitude control 
requirements must be provided by the propulsion system in the low »• 
speed flight regime and by both the propulsion system and the aircraft 
aerodynamic control surfacei? in the transitional flight regime until 
sufficient flight speed hac been achieved for full aerodynamic control. 
Integrated aircraft-engine controls caa be expected to be an essential 
requirement for achieving effective aircraft stability and control in 
the low flight speed and transitional flight regimes. 

The current V/STOL Propulsion Control Analysis Program represents 
Phase II of a multi-phase program for the development of the propulsion 
control technologies for achieving integrated aircraft-engine controls. 

The overall program is aimed at the development of the propulsion control 
design technology base, design analysis procedures, control logic require- 
ments, .and an overall system evaluation capability through piloted aircraft 
engine simulations. The initial Phase I program involved the definition 
of typical propulsion control requirements, the establishment of a long- 
range technology development plan, and the development of typical propul- 
sion control logic requirements for the vertical flight operating regime. 
The results of the Phase I program were published in Reference 1. The 
current Phase II program represents an extension of the Phase I effort 
to the transitional flight regime between vertical and horizontal flight. 
Subsequent program phases would further extend the program to the 
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conventional flight regime and would addreaa the overall problems of 
aircraft-engine control Integration. 

« 

The General Electric RALS (Remote Augmentor Lift System) concept 
was selected as baseline engine in the Phase I study In order to provide 
information on typical V/STOL propulsion system operational character- 
istics and control requirements. The baseline RALS engine was retained 
for the current Phase 11 studies so that the Phase X results could be 
used to establish Interface requirements with the vertical flight regime. 

Specific operational requirements were established for the following 
transitional processes: 

• Vertical takeoff at maximum weight to horizontal accelerated 
climb at maximum power setting.. 

m 

• Horizontal descent at flight idle to vertical landing at 
minimum weight. 

A number of specific operating points were identified along each 
of the above trajectories and corresponding steady-state operating 
requirements were established. Control mode studies were, conducted at 
each individual operating point to define closed>*loop control require- 
ments. Linear state-space models were developed' for each point and 
were used to define a multi-variable regulator design using the K-Q 
Matrix technique. Regulator gain schedules were obtained by regression 
of the individual designs. One set of schedules was obtained for the 
vertical-to-horizontal transition and a second set for the horizontal- 
to-vertlcal transition. 

A non-linear transient model of the baseline engine was constructed 
from regressions of the individual engine component operating character- 
istics. The component regressions were combined with regressions of 
variable gas property thermodynamics and with the non-linear models of 
the augmentors, nozzles, sensors, and actuators developed under the 
Phase I program to obtain a non-linear transient model of the propulsion 
system. The propulsion system model was used in conjunction with a 
mathematical model of the control logic to conduct simulation studies 


of typical transition operations. Additional mathematical models of 
inlet and nozzle effects and of a typical RALS aircraft configuration 
were also developed. The inlet and nozzle effects model defines the 
effects of high inlet distortion and nozzle deflections on the engine 
and will be combined with the propulsion system model during a subse~ 
quent phase of this program. The aircraft simulation has been used 
to determine propulsion control requirements and operating charactes" 
Istics during the transition process. 

This report contains a description of typical steady-state operat- 
ing requirements during the transition process » a description of 
recommended closed-loop control configurations and corresponding regu- 
lator gain requirements ,1 and a discussion of the component regression 
process for the development of the engine model and a brief description 
of the resulting model. It also contains an overall propulsion control 
concept for integrating engine, flight, and transition control require- 
ments and a revised long-range technology plan for achieving a piloted 
simulation capability, 

3.0 BASELINE PROPULSION SYSTEM AND CONTROL REQUIREMENTS . 

3.1 VCE/RALS Description 

V/STOL propulsion systems are required to provide thrust, bleed air, 
and power extraction in each of the following operational modes: 

• Vertical Mode - Requires vertical lift thrust over the range 
of takeoff and landing weights to achieve ascent, descent, and 
hover. 

• Attitude Control - Requires direct thrust or bleed air for 
achieving roll, pitch, and yaw control at low speeds when 
aircraft aerodynamic control surfaces are ineffective, 

• Transition Mode - Requires a continuous thrust vector rotation 
capability for achieving transitions between vertical and con- 
ventional horizontal flight. 

• Horizontal Flight Mode - Requires propulsive thrust, customer 
bleed, and power extraction for conventional horizontal flight. 

In order to achieve a high performance V/STOL capability, it will he 
essential to integrate each of the above operating mode requirements 
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into a vuciabXe cycle engine concept which can be adapted to each set of 
requirements and which minimizes the need for special geometry for any 
one mode of operation* 

The General Electric Variable Cycle Enginc/Remote Augmentor Lift 
System (VCE/HALS) was conceived for this purpose and was selected as the 

* baseline engine for the previous Phase I Study of V/STOL Propulsioi'' 
Control. It has been retained as the baseline propulsion system for £he 

* current Phase II study and is illustrated in Figure 1. The VCE consists 
of a single-bypass turbofan engine with the following features: 

• The fan is split into a forward two-stage block driven by the 
livw pressure turbine and a rear single-stage block driven by 
the high pressure turbine. 

• Variable rear block fan. stators and a variable low pressure 

turbine nozzle (VALPTN) for internal flow and pressurj? ratio 
control. * 

• A variable area bypass injector (VABI) for mixing fan and 
turbine discharge flows into the mixed-flow augmentor. 

• An Augmented Deflector Exhaust Nozzle (ADEN) for varying 
the primary thrust vector angle. 

The EALS concept adds to this VCE a remote augmentor and a vectorable 
remote nozzle which can be located near the nose of the aircraft in 
order to maximize the separation distance between the two thrust vectors. 
The remote system is operated by extracting most of the bypass air through 
a manifold and ducting it forward to the remote augmentor and nozzle. 

The remote system is shut down by modulating the VABI to direct all of 
the fan discharge air to the mixed-flow augmentor and ADEN nozzle. 

The current VCE/RALS concept uses both augmentors and nozzles in con- 
junction with engine speed modulation for vertical thrust control, 
mixed-flow and remote augmentor fuel-flow modulation for attitude con- 
trol, remote flow modulation for transition trajectory control, and 
only the basic VCE for horizontal thrust control. 

3.2 Propulsion Control For Vertical Operation 

The propulsion control concept for the vertical operating mode was 
developed in the Phase I study of the current program and is illustrated 
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Figure 1. VCE/RALS Baseline Engine 



in Figure 2. It consists of a multi-variable regulator for controlling 
the total engine thrust and a set of feed-forward schedules for modu- 
lating augmentor fuel-flows in response to attitude control demands. 

The multi-variable regulator operates over the range of 60-100% 
takeoff thrust to provide total thrust control during vertical takeoff 
and climb, hover, and vertical descent and landing. It consists of 
four closed control loops which set primary fuel flow, low pressure 
turbine nozzle area, and primary and remote Oixhaust nozzle areas to 
achieve zero feedback errors in the fan speed, turbine temperature, 
compressor discharge pressure and fan discharge duct Mach number loops. 
Nominal mixed-flow and remote augmentor fuel-flows are scheduled open- 
loop along with the remaining VCE geometries (rear-block faiL stators 
and VABI area). The regulator responds to a total thrust demand signal 
from the aircraft flight computer and represents the primary fuel con- 
trol system. It must be Integrated with accel/decel fuel schedules, 
limit protection schedules, bleed and power extraction compensation, 
etc. and these will be addressed during a subsequent phase of this 
program. 

4. 

The feed-forward schedules operate over the range of - 12% of 
nominal ADEN and remote nozzle thrust and provide the response to atti- 
tude control demands from the aircraft flight computer. These schedules 
provide mixed-flow and remote augmentor fuel flow modulation and ADEN 
and remote nozzle thrust deflections. The current VCE/RALS concept 
provides total thrust modulation for height control, differential 
thrust modulation for pitch control, differential thrust deflection 
for yaw control, and customer bleed extraction for roll control. It 
has been assumed that the aircraft flight computer demand signals will 
be the actual thrust magnitude and direction requirements. 

The thrust and attitude control systems interact through the aug- 
mentor fuel flows and exhaust nozzle areas which are common to both 
systems. The interaction has, however, been minimized by the addition 
of gain compensation schedules to both the regulator and the feed-forward 
schedules. Regulator gain compensation modulates augmentor fuel flow 



Figure 2. Propulsion Control for Vertical Operations. 










and nozzle areas at constant thrust to provide gas generator regu- 
lation with no effect on attitude control. Conversely » feed forward 
compensation modulates augmentor fuel flow and nozzle areas Co provide 
nozzle thrust regulation with no effect on gas generator speeds or 
temperatures. 


3 , 3 Flight Transition Control Requirements 

Two types of flight transitions between vertical and horizontal 
flight liave been examined: 

• VerCical-to-horizontal (VTH) transition from Vertical Takeoff 
to horizontal Climb rating. 

• Horizontal-to-Vertlcal (HTV) transition from Flight Idle Descent 
(FID) to Vertical Landing rating, 

The VTH transition involves the rotation of the ADEN and remote nozzle 
thrust vectors towards the horizontal, flow-split modulation to reduce 
the flow to the remote system, and the shut-down of Che remote augmentor. 
The HTV transition Involves the above events in reverse order plus a 
throttle burst from FID to Landing and light-offs of both the mixed-flow 
and remote augmentorsi Both transitions require the propulsion system 
to meet aircraft attitude control requirements during the low speed 
flight regime. * 

The following propulsion control requirements have been established 
for the flight transition phase; 

• Filot control of Power Lever Angle (PLA) and Thrust Vector Angle 
(TVA) in order to control the shape of the transition trajectory 
and its elapsed time. 

• Aircraft flight control computer to provide individual PLA 
and TVA demands to the propulsion control system. 

• Attitude control demands would be in the form of thrust magni- 
tude and vector angle corrections to both the /DEN and remote 
nozzles. 

• TVA demand establishes the applicable propulsion control regime - 
horizontal control for 0“, vertical control for 90“ TVA, and 
flight transition control for all intermediate values, 

• Stable mixed-flovj and remote augmentor operation must be available 
over a broad range of PLA and TVA demands. Augmentor Start-up 


and shut-down must be scheduled as a function of TVA (and 
possibly PLA as well). 

• Flow spilt between the VABI and the remote system must be 
scheduled as a function of TVA to minimize pitching moments 
due to the propulsion system. 

• Nominal mixed^flow and remote augmentor fuel schedules must 
be used during flight transition in order to reserve the full 
- 12% thrust modulation capability for potential attitude 
control needs. 

The above requirements have been used as the basis for the flight tran- 
sition control concept described in the subsequent sections of this 
report. Sea level static steady-state operating conditions were calcu- 
lated for a number of operating points along a typical VTH transition 
and a corresponding HTV transition. These data are described in Section 
4. These data were then used as the basis for the propulsion control 
regulator designs described in Section 5. 

4.0 STEADY-STATE OPERATION 


4.1 Horizontal Operating Mode 

Steady-state studies of the horizontal operating mode were .conducted 
in order to establish engine operating requirements for achieving minimum 
SFC over the full-throttle range from Idle to Intermediate Power Setting. 
Table 1 summarizes the results achieved for operation at SLS/Std. + 31®P 
and illustrates the variation in control parameters, potential sensed 
parameters, and performance parameters. These data were used to establish 
the starting point for the Horizontal-to-Vertical (HTV) transition and 
the end point for the corresponding Vertical-to-Horizontal (VTH) transition. 

4.2 Vertical-to-Horizontal Transition 

Figure 3 illustrates a typical VTH transition which contains the 
following key operating points: 

1) Hover or vertical climb with the remote flow control valve 
(RFCV) full open, 100% of the bypass duct flow going to the 
RALS nozzle, and the RALS system providing 45% of the total 
nominal thrust. The ADEN nozzle is in the vertical position 
(90°). 
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Figure 3. Typical Vertical-To-Horizontal Flight Transition. 



2) Start of transition v/lth the ADEN nozzle rotated to a 
thrust vector angle (TVA) of 77.4“ from horizontal. 

Remote flow transfer has been initiated by opening the 
reer VABI and the RALS thrust has been reduced to 
maintain a zero propulsive pitching moment on the 
aircraft. The RALS system is operating at 99% of 
nominal flow and producing 44% of the total thrust. 

3) TVA rotation is decreased to 35.7“ and the rear VABI 
is opened to reduce RALS flow to 80% nominal. The 
RFCV must be partially closed to maintain RALS burner 
Mach number. RALS thrust has been reduced to 32% and 
air speed and aerodynamic lift are increasing. 

4) The TVA has been reduced to 10.8“ and the remote flow 
to 50% of nominal. Increased throttling loss through 
the RFGV and reduced augmentor fuel flows reduce the 
RALS thrust to 13% of nominal. Air speed and aerodynamic 
lift continues to increase. 

5) TVA has been reduced to 2.3“ and the remote flow to 30% 

of nominal. This is the minimum flow rate at which the 

RALS burner can be maintained and the burner AT has been 
reduced to a minimum of 150“F in preparation for shut-off. 

6) TVA has been reduced to 2.1“ and the remote burner has 

been shut-off. 

7) The RFCV has been completely shut-off and all duct flow 
passes through the rear VABI and ADEN nozzle completing 
the transition. 

Table 2 contains a summary of the control, manipulated and performance 
variables corresponding to the above points along the VTH transition. 
Note that there is relatively little change in engine speed throughout 
the VTH transition and that the primary augmentor has been left on in 
order to achieve a fast acceleration to cruise speed and altitude. 


4.3 Horizontal-to-Vertlcal Transition 

The HTV transition is expected to be the most critical phase of the 
VCE/RALS operating envelope in that it requires a transfer of about 60% 
of the bypass duct flow to the RALS system before igniting the RALS 
burner and a substantial throttle transient from Flight Idle to at least 
60% of Nominal Vertical Takeoff Thrust. 


Figure 4 illustrates a typical HTV transition containing the follow- 
ing operating points: 
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TABLE 2. TYPICAL VERTICAL-TO-HORIZONTAL OPERATING CONDITIONS 
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Figure 4. Typical Horizontal ~To-Ve;rtical Flight Transition. 



! 

1) Flight Idle Descent corresponding to low power horizontal 
operation with no HALS flow and both augmentors turned off, 

The RFCV is closed and the ADEN TVA is at 0* . 

i 

2) TVA has been increased to 3.9" and the RFCV has been opened 
to transfer 30^ of the duct flow to the RALS' system. Both 
augmentors are turned off and the RALS thrust is 5% of total 
thrust. 

3) TVA has been increased to 11.5" and the RALS flow has been 
Increased to 60% of the total flow. The RALS burner Mach 
number is adequate for burner ignition. RALS thrust is 14% 
of the total thrust. 

4 & 5) The primary ai,d RALS augmentors are ignited consecutively 
and brought to minimum augmentation at 150"F temperature 
rise. 

6) The ADEN TVA has been increased to 41.2" and the RALS flow 
has been increased to 80% of the total flow. Fan speed has 
been Increased and both augmentor fuel flows have been raised 
to achieve a RALS thrust of 35% of total. The RFCV has been 
opened in conjunction with the RALS nozzle area to maintain 
the required RALS burner Mach number. 

7) The RFCV is fully open and the rear VABI is closed. RALS 
flow is 100% of the duct flow and the ADEN TVA is at 90" 
from horizontal. Fan speed has been increased to achieve 
60% of the nominal VTO thrust and the ^LS thrust is 45% 

of the total completing the transition into hover or vertical 
landing. 

Table 3 contains a summary of the control, manipulated, and performance 
variable variations corresponding to the above HTV transition. 



Figure 5 illustrates a typical flight transition regime at sea level, 

ADEN gross thrust has been plotted against the ADEN Thrust Vector Angle. 

Horizontal operation is shown from Flight Idle to Maximum Augmentation at 
0* TVA. Vertical operation is shown from Minimum Weight Hover (60% 

Nominal Takeoff) to Maximum Weight Hover (100% Nominal Takeoff) at 90" | 

TVA. The upper bound represents the nominal VTH transition described in - | 

Section 4.3 and the lower bound the nominal HTV transition. Lines of I 

constant RALS thrust are Indicated as functions of the ADEN thrust and | 

TVA and correspond to zero propulsive pitch operation (no pitch rate i 

I 

demand). RALS and ADEN thrust modulation and deflection are available | 

I 

for attitude control by modulating the primary and RALS augmentor fuel flows. * 
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TABLE 3. TYPICAL HORIZONTAL-TO-VERTICAL OPERATING CONDITIONS 
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The transition process and trajectory shape are controlled by the 
individual pilot demands: 

• Power Lever Angle which sets the engine fan speed and the total 
thrust at nominal augmentation levels. 

• Thrust Vector Angle which sets the nominal ADEN deflection and 
the nominal RALS thrust to achieve a zero propulsive pitch rate. 

TVA determines the relative accelerations parallel and petpendi- 
cular to the aircraft and^ therefore, the rate of rotation of 
the velocity vector. 

• Pitch Rate which modulates the nominal RALS thrust for zero pro- 
pulsive pitch in order to achieve the required pitch rate. The 
Pitch Race demand, therefore » controls angle of attack during the 
transition process. 

Note that the nominal RALS and i/^EN thrust requirements can be deter- 
mined uniquely from the total Gross Thrust and Thrust Vector Angle demands. 
The thrust modulation requirement is defined by the Pitch Rate demand. 

Thrust modulation capabilities available for height and attitude 

control are illustrated in Figure 6. The nominal thrust modulation caps- 
"f* 

bility is -.12% for both the RALS and ADEN nozzles throughout most of the 
flight transition regime. At maximum weight takeoff, the minimum ADEN 
temperature is only 765“R below the nominal temperature limiting the low 
side modulation to 8.9%. This capability could^, however, be increased 
by utilizing internal speed or geometry variations in this part of the 
operating regime, if necessary. 

4.5 Remote Flow Control Valve 

The RALS burner is required to operate over a broad range of flow 
rates, pressures, and temperatures throughout the flight transition regime. 
Stable operation of the RALS burner requires an inlet Mach number in the 
range of .13 - .15 and the RALS geometry has been sized to operate in this 
range at nominal remote flow rates. The rear VABI controls the flow split 
between the RALS and ADEN systems and the RALS nozzle area provides the 
means for controlling the burner Mach number. 

Operation of the Remote Flow Control Valve is required in the low 
flow regime associated with the mid-transition region in order to throttle 
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the fan discharge pressure to maintain an acceptable burner Inlet Mach 
number. Figure 7 Illustrates the effect of the RFCV during the VTK 
transition of Table 2. Typical operation involves the following operating 
points: 


1) At nominal remote flow rate, the RFCV is open and the Mach 
number is controlled by the RALS nozzle area. The pressure 
drop between fan discharge and the burner inlet is a minimum 
and is due to friction and the head loss associated with 
turning the flow forward. 

2) At 80% of nominal flow rate, both the RFCV and RALS areas are 
reduced. Burner inlet Mach number is maintained by the addi- 
tional throttling losses in the system. Mote that a large 
reduction in RFCV area is required in this regime as the 
throttling losses are relatively insensitive to the area change. 

3) The RFCV is choked at about 50% of nominal flow and assumes the 
flow control function from the rear VABI. The RALS nozzle area 
determines the dump pressure behind the RFCV and therefore its 
losses. RALS burner pressure is substantially more sensitive to 
the RFCV and mass flow rate at this point. 

4) This point represents the minimum mass flow rate at which the 
required burner inlet Mach number can be maintained due to 
external ambient pressure limitations. The step change in 

. RALS nozzle area corresponds to RALS burner shut-down from 
operation at minimum temperature rise (150“F AT). 

5) The RFCV is closed and the duct flow hqs been completely transi- 
tioned from the RALS system to the ADEN nozzle. 

4.6 Component Operating Characteristics 

* 

Component operating characteristics from the horizontal operating 
points of Table 1 and the flight transition points of Tables 2 and 3 have 
been used in the development of the simplified component-level transient 
model of the baseline engine described in Section 6. These data have been 
supplemented by additional steady-state operating data involving off-nominal 
operating schedules. Figures 8 and 9 illustrate the front-block fan and 
rear-block fan operating pressure ratio variations with corrected flow 
corresponding to the above operating data. The associated rear-block 
fan stator variation is also illustrated in Figure 9. Figure 10 contains 
similar data showing the compressor pressure ratio variation with corrected 
flow. Figures 11 and 12 contain the corresponding high pressure and low 
pressure turbine energy variations with turbine inlet corrected speed. Low 
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pressure turbine variable stator position variation Is also* shown In 
Figure 12. 


5.0 PROPULSION CONTROL LOGIC 


5.1 Control Mode Studies 

Control mode studies were performed for each of the VTH and HIV 
operating points In Tables 2 and 3 In order to establish closed-’loop 
control configurations for the subsequent regulator design activity. 

The mode studies were based upon the following: 

• Component quality, control tolerance, and deterioration 
models described In Reference 1. 

• Partial derivative models for each of the Individual VTH 
and HTV operating points derived from the baseline engine 
steady-state cycle deck. 

• The selected control mode configuration for the vertical 
operating regime at Takeoff (100% Nominal Thrust) and 
Landing (60% Nominal Thrust). 

• Additional candidate control modes where appropriate. 

•V . * 

Table 4 summarizes the results of the VTH mode studies. These 
results Indicated the following: 

• A four closed-loop control configuration should be retained over 
the total VTH transition process. 

• Closed-loop control (sensed)^ variables should be Fan Speed 
(PCN2) , Low Pressure Turbine Discharge Temperature (T5) , 
Compressor Discharge Static Pressure (PS3), and Bypass Duct 
Mach Number (XM93). 

• Closed-loop manipulated variables over the total VTH transition 
should be ADEN Nozzle Area (Ab) , Primary Fuel Flow (WF36), and 
Low Pressure Turbine Variable Stator Position (STP49). The 
fourth closed-loop manipulated variable should be RALS Nozzle 
Area (A88) In the high remote flow regime and Rear-Block Fan 
Variable Stator Position (STP22) In the low remote flow regime. 

• Open-Loop manipulated variables are the rear VABI area (A27), 
nominal Primary Augmentor Fuel Flow (WF6), and nominal RALS 
Fuel Flow (WF86) . STP22 would be controlled open-loop at high 
flow rates and A88 at low flow rates. 
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Control of the RFGV in the mid-transition region was assumed to be 
open-loop but: this assumption should be reviewed in more detail during 
subsequent studies* 

Table 4 contains the expected variations in thrust » stall margin, 
temperaCure, and fuel flow corresponding Co the selected control modes. 
These data contain no significant differences from the prior mode 
studies for the vertical operating mode. 

Additional mode studies Were conducted for 3-loop and 5-loop 
systems. The best 3-loop modes Involved unacceptable performance 
variations and the best 5-loop modes did not provide sufficient improve- 
ments to warrant the added complexity of a fifth loop. 


Similar mode studies Were conducted for the HTV transition points 
of Table 3 and the results are summarized in Table 5. These results 
indicated the following: 


• A single closed-loop control was acceptable for operation at 
Flight Idle Descent (FID) . This control should use PCN2 as 
the sensed variable and WF36 as the closed-loop manipulated 
variable. All other manipulated variables should be scheduled 
open-loop . 

• A four closed-loop configuration was required to achieve 
acceptable performance variations for each of the remaining 
operating points along the HTV transition. The HTV transition 
could use the same closed-loop control and manipulated variables 
as the VTH transition process. 


5.2 Controller Design 

The Controller design process followed the same procedure used in 
the prior Phase I program (Reference 1) and involved the following: 

• A rectangular matrix of balanced steady-state partial derivatives 
was used to establish a state-space model with two states, four 
inputs, and four outputs for each of the VTH and HTV operating 
points. The state-space models were scaled and structured in the 
ABCD format. Figure 13 contains a schematic of the Linear Engine 
state-space model format. 
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TABLE 4 - VERTICAL TO HORIZONTAL TRANSITION MODE STUDY 


OPERATING POINT 

T/0 

S/T 

A/RFCV 

■OBH 

M/R 

IBHI 

CLOSED 

PCN2 

(A8) 

(A8) 

(A8) 

(A8) 

(A8) 

(A8) 

LOOP 

T5 

(WF36) 

(WF36) 

(WF36) 

(WF36) 

(WF36) 

(WF36) 

CONTROL 

PS3 

(STP49) 

(STP49) 

(STP49) 

(STP49) 

(STP49) 

(STP49) 

VARIABLES 

XM93 

(A88)* 

(A88)*_. 

JSTP22)* 

(STP22)* 

(STP22)'* 

(STP22)* 

CLOSED 

A8 

(PCN2) 

(PCN2) 

(PCN2) 

(PCN2) 

(PCN2) 

(PCN2) 

LOOP 

WF36 

(T5) 

(T5) 

(T5) 

(T5) 

(T5) 

(T5) 

MANIPULATED 

STP49 

(PS 3) 

(PS3) 

(PS3) 

(PS3) 

(PS3) 

(PS3) 

VARIABLES 

STP22 

A88 

(XM93)** 

y 

(XM93)*^ 

^(XM93) 

(XM93) 

(XM93) 

(XM93) 

MINIMUM 

|AF689 

-3.28 

-3.74 

-5.66 

-9.4 

-7,91 



?AFG9 

-6.87 

-6.67 

-4.93 

-3.37 

-2.17 

-4.28 


^ASM2 

-1.79 

-1.79 

-1.79 

-1 .79 

-1.79 

-1.84 


“^SM22 

-4.84 

-4.81 

-3.04 . 

-0.19 

-3.55 

-0.11 


%AsM25 

-2.14 

-2.14 * 

-2.13 

-2.15 

-7.10 

-2.22 

MAXIMUM 

AT41(“For “R) 

m 

47 

48 . ' 

46 

36 

43 


AT48(°For °R) 


40 

40 

40 

41 

38 


Asfc 


.051 

' .044 

.042 

.024 

.10 


AWF36 (PPH) 

Bi 

240 

218 

241 

215 

229 


T/O-Takeoff 
S/T-Start Transition 

A/RFCV-Activate Remote Flow Control Valve 
M/T-Mid Transition 
M/R-Min RALS (AT = 150'’) 

MH -Max Horizontal 

* -Manipulated Variable Associated With XM93 For Operating Point 
** -Control Variable Associated With A88 For Operating Point 
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TABLE 5 - HORIZONTAL -TO-VERT I CAL TRANSITION MODE STUDY 


OPERATING POINT 

FID 

30 RF 

60 RF 

RALS 

80 RF 

L 

CLOSED 

PCN2 

(WF36)* 

(A8)* 

m* 

(A8)* 

(A8)* 

(A8)* 

LOOP 

T5 


•KWF36) 

(WF36) 

(WF36) 

(WF36) 

(WF36) 

CONTROL 

PS3 


(STP49) 

(STP49) 

(STP49) 

(STP49) 

(STP49) 

VARI/iBLES 

XM93 


{STP22) 

(STP22) 

(STP22)- 

it<A88) 

(A88) 

CLOSED 

A8 


.-WPCN2) 

(PCN2) 

(PCN2) 

(PCN2) 

(PCN2) 

LOOP 

WF36 

(PCN2)** 

(T5)** 

(T5)** 

(T5)** 

(T5)** 

(T5)** 

MANIPULATED 

STP49 


(PS3) 

(PS3) 

(PS3) 

(PS3) 

(PS3) 

VARIABLES 

STP22 


(XM93) 

(XM93) 

(XM93)^ 




A88 




> 

^XM93) 

(XM93) 

MINIMUM 

%AfG89 

immum 

-6.88 

-7.88 

-7.53 

-4.00 

-4.28 


%AfG9 

-3.5 

-2.42 

-L58 

-2.63 

-5.06 

-6.36 


%AsM2 

-1.9 

-1.48 

-1.5 

-1.5 

-1.76 

-1.84 


%AsM22 

-2.99 

-1.94 

-2.71 

-2.55 

-3.19 

-4.00 


%AsM25 

^4.17’ 

-10.19 

-10. 

-9.98 

-9.81 

-7.09 

MAXIMUM 

AT41(°For °R) 

82 

22 

f 

23 

23 

31 

29 


AT48(°Ror °R) 

76 

23 

24 

25 

29 

33 


Asfc 

.06 

.03 

,035 

.043 

.07 

.055 


AwF36 (PPH) 

95 

33 

39 

40 

82 

79 


FID - Flight Idle Descent 
30 RF - 30% RALS Flow 
60 RF - 60% RALS Flow 
RALS - RALS (AT = 150°) 

80 RF - 80% RALS Flow 
L - Landing (60% Max. T/0 Power) 

* - Manipulated Variable Associated With PCN2 For Operating Point 
**-Control Variable Associated With WF36 For Operating Point 
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AXN2 

AT5 

APS3 

AXM93 


Figure 13. Four Imput-Four Output Linear Model (LEM) 
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• The K/Q tMCKlK technique was then used to design Proportional/ 
Integral Controllers (PIC) for each individual operating point. 

The design technique was applied to the open-loop structure 
illustrated in Figure 14. Identical scale factors were used 
(or all operating points. 

• A computer generated closed-loop diagram of the system was 
obtained and is illustrated in Figure 15. Note that the 
diagram indlcotes the inter-relationships between the Linear 
Engine Model (SYS), Actuator Matrix (PREl), Controller Matrix 
(PUE2) , and the Sensor Matrix (F/Bl) , Figure 16 contains a 
more conventional diagram of the closed-loop structure in the 
standard multi-variable form. 

The resulting regulator comprised the following elements: 

• First Precompensator (PREl) - A 4x4 diagonal matrix containing 
the dynamics of the actuators. PREl was the same for all 
transition operating points. 

• Feedback Compensators (F/Bl) - A 4x4 diagonal matrix containing 
the dynamics of the sensors. F/Bl was the same for all transition 
operating points. 

• Controller and Second Precompensator (PRE2) - A 4x4 matrix con- 
taining PrOportional/Integral Controller (PIG) Laplace Transfer 
functions for each of its 16 elements. Different transfer 
functions were obtained for each individual operating point in 
the VTH and UTV transitions. Table 6 contains the VTH results 
and Table 7 the results for the HTV transition. Note that Flight 
Idle Descent utilizes a single-loop control and, consequently, 
its PRE2 matrix contains only a single set of transfer functions, 

✓ 

Appendix A contains the results of a series of time domain plots for 
unit step demands applied simultaneously to all Inputs of the closed-loop 
regulator for each of the VTH and HTV operating conditions. In all cases 
ctoss-coupllng occurred within the first quartet Second and steady-state 
was achieved within two seconds. 


5.3 Gain Schedules 

The proportional and integral gain constants from Tables 6 and 7 
were then represented by linear regression as a function of an arbitrary 
flight control demand (FCD) parameter. F'CD values were determined in 
order to achieve the best fit for the major diagonal coefficients. 
Although FCD is a function of TVA and Power Lever Angle, no attempt has 
been made to define the functional relationship. Linear regressions were 
obtained for each of the 32 coefficients for the VTH transition and for 
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PROPORTIONAL AND INTEGRAL GAIN CONSTANTS FOR THE 
VERTICAL-TO-HORIZOMTAL. TRANSITION CONTROLLER 


INTEGRAL 

GAIN 

COEFFICIENTS 


7Y 

POINT 

TAKEOFF 

TK02F 

START 
TRANS ITIOI 
VH22P 

ACTIVATE 

RFCV 

VH32P 

HID 

TPANSITIOI 

VH42P 

HIN RALS 
VH52P 

MAX 

HORIZONTA 

MXH2P 

'ED VARIABLE 

A88 

A88 

STP22 

STP22 : 

STP22 

STP22 

all 

1.51 

1.42 

0.92 

1.44 ■ 

1,53 

.13 

a12 

-0.69 

-0.63 

-0.53 

-0.49 

-0.59 

-0.94 

a13 

0.91 

0.86 

1.67 

0.91 

! 0.94 

1.21 

al4 

0.59 . 

0.54 

-0.01 

0.05 

0.07 

-0.08 

a21 

0.64 

■ 0.63 

0.42 

0.64 

0.65 

0.12 

a22 

0.49 

0.48 

0.53 

0.59 

0.63 

0.62 

a23 

0.39 

0.39 

0.68 

0.40 

0.38 

0.35 

a24 

0.21 

0.20 

-0.01 

0.005 

0.008 

-0.03 

a31 

-0.43 

-0.42 

-0.86 

-0.36 

-0.35 

-0.56 

a32 

-0.47 

-0.46 

-0.36 

-0.30 

-0.32 

-0.34 

a33 

0.96 ■ 

0.96 

1.56 

0.87 

0.82 

0,79 

a34 

0.59 

0.58 

0.06 

0.11- 

0.12 

0.03 

. a41 

-0.09 

-0.11 

1,30 

3.20 

3.11 “ 

2.13 

1 

a42 

0.34 

0.35 

-1.55 

-1.98 ■ 

-2.09 

-2.36 

a43 

-0.19 

-0.20 

1.24 

2.17 

2.11 

2.30 

a44 

0.33 

0.33 

-0.63 . 

^0.53 

-0.50 

-0.86 

bn 

7.08 

6.58 

8.43 

7.86 

8.62 

17.55 

b12 

-0.62 

-0.57 

'-0.37 

-0.32 

-0.39 

-0.68 

bl3 i 

-0.54 

-0.51'' 

-3.42 

-3.27 

-3.37 

-4.92 

bU 

-0.54 

-0.40 

0.85 

0.80 

0.94 

1.3c 

b2I 

1.79 

1.79 

2.07 

1.35 

1.33 

2.67 

bZZ 

0,61 

0.60 

0,73 

0.80 

0.83 

0.84 

bZ'^ 

3.99 

3.99 

3.21 

3.65 

3.50 

3.39 

b24 

-0.39 

-0.38 

0.14 

0.08 

0.07 

0.20 

b31 

-1.88 

-1.85 

-1.09 

-2.51 

-2.77 

-6.76 

b32 

-0.37 

-0,36 

-0.17 ! 

-0.13 

-0.15 

-0.17 

b33 

3.12 

3.07 

1.22 

2.02 

2.14 

2.51 

b34 

1.36 

1.32 

2.34 

2.12 

2,20 ; 

2.16 

b41 

4,94 

4.86 

-8.42 

7.75 

8.66 

8.95 

b42 

0.32 

0.32 

-1,37 

-1.65 

-1.75 

-1.98 

b43 

-7.89 

-7.86 

18.95 

7.27 

6.90 

5.34 

b44 

4.45 

4.47 

-12.85 

-11.50 

-11.24 

-13.43 


OPERATING POINT 


SPECIAL MANIPULATED VARIA3LE 


PROPOTlTIOfWL 

GAIN 

COEFFICIENTS 


INTEGRAL 

GAIN 

COEFFICIENTS 


STP22 

STP22 

STP22 

ASS 

1.08 

0.76 

0.86 

0.90 

-0.77 

-0.58 

-0.61 

-1.10 

7.88 

5.08 

5.22 

3.84 

-0.0007 

-0.02 

-0.02 

0.61 

0.10 

0.09 1 

0.09 

0.14 

0.18 

0.19 

0.19 

0.19 

0.36 

0.36 

0.35 

0.54 

-0.005 

-0.007 

-0.006 

0.09 

-14.06 

-9.59 

-9.56 

-4.43 

-0.43 

-0.47 

-0.46 

-1.16 

6.91 

5.30 

5.18 

4.60 

0.20 

0.14 

0.14 

0.75 

2.97 

2.67 

2.90 

-1.20 

-4.19 

-3.88 

-3.83 

1.20 

26.81 

20.43 

20.38' 

-1 .21 

-1.27 

-1.23 

-1 .20 

0.40 

5.27 

4.86 . 

' 5.03 

5.92 

-0.26 

-0.22 

-0.24 

-0.89 

-7.12 

, -4.04 

-4.42 

0.33 

3.16, ■ 

2.16 

2.21 

-0.96 

0.39 

0.51 

0.50 

0.97 

0.25 

0.28 

0.28 

0.33 

2.15 

2.27 

2.30 

2.86 

0.10 

0.11 

0.11 

-0.07 

^34.74 

-28.44 

-28.51 

-12.27 

• 

0 

1 

-0.12 

-0.12 

-0.69 

21.26 

15.77 

15.86 

10.95 

5.44 

5.20 

5.16 

1.29 

17.94 

18.80 

19.70 

11.11 

-3.04 

-2.77 

-2.74 

0.91 

16.88 

13.57 

12.51 

-18.85 

-16.26 

-15.68 

-15.51 

9.51 


A88 


.07 

-2.30 

4.52 

1.81 

0.07 

0.32 

0.46 

0.17 

-3.68 

-1.56 

3.80 

1.49 

-1.04 

0.37 

-0.09 

0.38 


w I 
0.16 

I. 77 
0.45 
3.07 

-0.16 

-18.89 

- 1.20 

11,56 

2.85 

II. 26 
0.34 

-8.99 

5.01 











Qiich oi: fclvo 32 coe|fiGiai\t{5 Sm the Ul'V ttiinsiCiou (but with thu giuglu- 
loop FID ease omittud) . Tablua 8 and 9 contain tha taswltlng gains fot 
each VTU and H'i'V opavatlng points obtainad ftom tha roaults of tha 
fagtasaion ptocaas. Gompafison batwean Tablaa 6 and 8 and batwean 7 Jind 
9 ptevida an indication of t)\a magnituda of any cosuitlng fitting actors ^ 

Appandix 8 contains tha rosuXts of a aacond sarias of tima domain 
results for unit atop demands applied simultanaously to all regulator 
ittputa using the gains from I'ables 8 and 9, Iha major effect of the 
regression errors was to increase the amount of interaction (high frequency) 
hetween control loops, the response continued to achieve steady~state 
although in several cases as a\nch as five seconds was required, 1‘hts was 
particularly true for the control variable T5. 

I'he determination of the robustness of a multi~variable control 
system is one of the most dynamic aspects of advanced control theory and 
better methods for determining robustness continue to surface. Robustness 
tests are now available for parameter variations, gain variations, and 
frequency dependent variations for specific blocks of a closed**loop 
multi-variable control system. Robustness studies were conducted on the 
transition regulator for both sets of gain coefficients. As would be 
expected, the linear regression model was not as robust as its parent 
model but its margins \Wre, however, adequate. 
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CONSTANTS FOR THE VERTICAI.-TO-HORIZONTAL TRANSITION CONTROLLER 


rCD VALUE 


OPERATING POINT 


SPECIAL MANIPULATED VARIABLE 


all 


GAIN 

COEFFICIENTS 


INTEGRAL 

GAIN 

COEFFICIENTS 



O o- cvj 
Lu cvi 

o o 

< fw: 


ABB 


1.59 

-0.61 

0.90 

0.59 

0.69 

0.48 

0.416 

0.21 

-0.436 

-0.467 

1.021 

0.6 

-0.148 

0.460 

-0.315 

0.379 


5.618 
-0,558 
-0.281 
-0.563 
1.677 
0.592 
3,997 
-0.404 
- 1,204 
-0;367 
2.993 
1 .34 
3.549 
0.413 
•10.768 
4.861 


o 

^>4 

1/1 



•-I Cu 


1/) O. CM 

a: 

scNjcr 
O < .Ui 

•m* 

^ es;. ZZ'^ 


s; H" > p» 



o o. 
r>i a. 

— fvj 

X c £ =: 
c o X X 


AS8 

STP22 

STP22 

1.52 

1.06 

1.02 : 

-0.62 

-0.65 

-0.66 

0.93 

1.13 

1.15 

0.52 

0.09 

0.06 

0.66 

0.47 

0.46 

0.^9 

• 0.58 

0.58 

0.42 

0.436 

0.437 

0.19 

,026 

.013 

-0,447 

-0.512 

-0.517 

-0.451 

-0.352 

-0.343 

1.016 

.986 

,984 

0.537 

0.159 

0.12S 

0.163 

2.031 

2.187 

0.160 

-1.64 

-1.79 

-0.037 

1.632 

0 

1 .772 

0.256 

-0.485 

-0.547 

6.28> 

10.301 

10.636 

-0.546 

-0.475 

-0.469 

-0.709 

-3.278 

-3.492 

-0.372 

Q.712 

0.867 

1.705 

1.873 

1.887 

0.618 

0.771 

0.784 

3.93 

3.526 

3.493 

-0.34 

0,044 

0.076 

-1.492 

-3.217 

-3.361 

-0.341 

-6.189 . 

! -0.176 

2.877 

2.183 

2.125 

1.443 

2.063 

2.115 

3.712 

4.687 

4.763 

0.157 

-1 .382 

-1.510 

-8.397 

5.829 

7.014 

2.793 

-9.615 

-10.649 




0.98 

- 0.66 

1.16 

0.02 

0.44 

0.59 

0.438 

-2.2x10' 

-0,523 

-0.335 

.981 

0.095 

2.342 

-1.94 

1.911 

-0.603 


10.97 
-0.463 
-3.707 
0.962 
1.901 
0.797 
3.459 
0.108 
-3.505 
-0.164 
2.067 
2.166 
4.850 
-1 .638 
8.2 

■11 .633 


2.308 

0.439 
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USE OF LINEAR REGRESSION TO DETERMINE PROPORTIONAL AND INTEGRAL GA I 


CONSTANTS FOR THE HORIZONTAL-TO-VERTICAL TRANSITION CONTROLLER 


FCD VALUE 


OPERATING POINT 



SPECIAL MANIPULATED VAP.IASLE 


STP22 

STP22 

STP22 

A88 


all 

3.76 

.98 

.92 

.86 

.48 

.42 

a12 


-0.51 

-0.65 

-0.79 

-1 .64 

-1.78 

a13 


6.407 

6.133 

5.858 

4.209 

3.934 

a14 


-0.187 

-0.021 

0.145 

1.138 

1.304 

a21 


0.095 

0.095 

0.096 

0.101 

0.102 

PROPORTIONAL ^22 


0.175 

0.184 

0.194 

0.253 

0.263 

GAIN ^23 


0.343 

0.361 

0.379 

0.485 

0.503 

COEFFICIENTS ^24 


-0.024 

-0.006 

0.012 

0.121 

0.139 

a31 


-12.150 

-11.179 

-10.207 

-4.378 

-3.406 

a32 


-0.337 

-0.456 

-0.576 

-1.295 

• 

cn 

a33 


6.101 

5.865 

5.629 

4.215 

3.979 

a34 


0.032 

0.161 

0.289 

1.060 

1.188 

a41 


3.299 

2.790 

.>.280 

-0.779, 

-1.2S9 

a42 


-4.507 

-3.896 

-3.286 

0.375 

0.935 

a43 


25.526 

22.460 

19.395 

1.002 

-2.063 

a44 


-1.422 

-1.212 

-1.003 

' 0.254 

0.464 

bll 

4.30 

5.067 

5.093 

5.118 

5.269 

5.294 

b12 


O 

O 

1 

-0.235 

-0,393 

-1 .343 

-1.501 

b13 


-6:337 

-5.290 

-4.243 

2.041 

3.083 

b14 


2.865 

2.485 

2.105 

-0.173 

-0.553 

b21 


0.334 

0.458 

0.581 

1.322 

1,445 

INTEGRAL 







b22 


0.251 

0.268 

0.285 

0.385 

■ 0.402 

GAIN 

b23 


2.143 

2.239 

2.336 

2.917 

3.014 

COEFFICIENTS 







b24 


0.134 

0.105 

0.076 

-0.09S 

-0.127 

b31 


-32.443 

-30.474 

-28.506 

-16,697 

-14.729 

b32 


-0.02 

-0.128 

-0.237 

-0.888 

-0.997 

b33 


18.650 

17.757 

16.865 

11.51 

10.618 

b34 


6.032 

5.571 

5.110 

2.344 

1.383 

b4I 


19.58 

18.626 

17.671 

11.544 

10.989 

b42 


-3.248 

-2.801 

-2.364 

0.328 

0.775 

b43 


17.515 

13.893 

10.270 

-11.467 

-15.090 

b44 


-18.368 

-15.423 

-12.477 

5.196 

8.142 





















5.4 Overall Control Concept 


Figure 17 contains a preliminary block diagram of the overall V/STOL 
Propulsion Control concept evolving for the baseline VCE/RALS propulsion 
system. It indicates control Inputs, outputs, and the principle inter- 
actions between the three major control sub-systems - the engine control, 
the thrust vector control and the flight control. 

The engine control provides the basic functions such as main fuel 
control. Internal variable geometry control, speed regulation, and engine 
limit protection. Engine control Inputs include external sensor data (P2, 
T2, Mach. No.). Power Lever Angle (PLA)demand, and closed-loop sensor 
data (PGN2, T5, PS3, and XM93) . The engine control logic includes: 

• Closed-loop demand schedules for PCN2, T5, PS3, and XM93 

• Open-loop demand schedules for STP22| and for STP49, A8 and A88 at 
reduced power operation 

• Gain schedules for establishing regulator proportional/integral 
gains as a function of Thrust Vector Angle (TVA) and PLA demands 
and possibly as a function of flight condition (P2, T2, Mach. No.) 

• Multivariable regulator Which converts the closed-loop error 
signals into demand signals for the closed-loop manipulated 
variables (WF36, STP49, A8, and A88) 

• Transition controls for limiting regulator outputs during large 
and/or fast throttle bursts and chops.' The transition controls 
Include fuel accel and decel schedules for stall and blowout 
protection; engine speed, temperature, and pressure protection; 
and any special logic required' for augmentor light-off and 
shut-down protection. 

The engine control logic provides output demands for primary fuel 
flow (WF36) and low pressure turbine stator position (STP49) to the 
position servos and nozzle area error signals (^A8R and Z1 a 88R) to the 
feed forward system. 

The thrust vector control sets the nominal thrust split between the 
RALS and ADEN nozzles and the ADEN thrust vector angle. Thrust vector 
control inputs include sensor data (T2, PCN2), TVA demand, and Pitch Rate 
(PR) demand. A Thrust Ratio (TR) demand could be used in place of indivi- 
dual TVA and PR demands and represents the ratio of RALS Gross Thrust to 
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Flight Control Thrust Vector Control Engine Control 



























ADEN Gross Thrust, The thrust vector control Includes: 

• Open-loop demand schedules for the rear VABI (A27) and the 
Remote Flow Control Valve (RFCV) which sets the flow split 
between the primary and remote augmentors. 

• The nominal ADEN deflectirn angle (ADL) demand which is set 
to the TVA demand. 

• RALS and ADEN nominal gross thrust demands (FG89D and FG9D) . 

ADEN gross thrust demand is dependent upon PLA demand, flow 
split demand, and the external operating condition (P2, T2, 

Hach. No.). RALS gross thrust demand is determined from the 
ADEN thrust demand and the thrust ratio demand implied by TVA 
and PR. Note that the individual thrust demands effectively 
set nominal augmentor temperatures in both the ADEN and RALS 
augmentors. 

The thrust vector control logic provides output demands for A27 and 
RFCV to the positon servos and nominal demands for ADL, FG89, and FG9 to 
the flight control system. 


The flight control responds to individual ADEN and RALS thrust magnl 
tude and deflection demands from the aircraft height and attitude control 
system. It includes the following elements: 

• ADEN deflection demand (^ADL) is added to the nominal demand 

signal from the thrust vector control (ADLD) to establish the 

actual demand. , * 

» 

• Thrust magnitude demands (AfG 89 and AfG 9) are added to the 
nominal demand signals from the 'thrust vector control to estab- 
lish the actual demands. 

• RALS deflection demands in the longitudinal direction (^RDL) and 
the transverse direction (ArDT) are used to set the actual RALS 
deflection demands. 

• The actual thrust demands and the nozzle area error signals from 
the engine control are input to the feed forward system to estab- 
lish augmentor fuel flow and nozzle area demands to the ADEN (WF6 
and A8) and the RALS (WF86 and A88) systems. Note that the feed 
forward system modifies all four demand signals to trim nozzle 
area errors at constant thrust and to trim thrust corrections at 
constant engine operating point. 

The flight Control logic provides demand signals directly to the corres- 
ponding positon servos. 


44 


Control logic has been defined under the Phase I and XI study 
efforts for the regulator gain schedules » the regulator* and the feed- 
forward system. Additional control logic must be developed under sub- 
sequent study programs to define the following: 

• Closed-loop and open-loop demand schedules 

• Thrust vector control demand schedules 

• Transition control schedules and logic 

• Gain schedules for the feed-forward system* if required 

Each of the above schedule requirements need to be examined over their 
respective operating regimes to determine scheduling requirements with PLA* 
TVA* and any external operating effects due to P2* T2* and flight Mach 
number. Control stability and robustness studies must be conducted at 
selected operating conditions to establish any additional control logic 
requirements for full range-full power operation. 

6.0 MATHEMATICAL MODELING 
6.1 Inlet and Nozzle Effects Model 

A realistic evaluation of V/STOL steady-state and transient performance 
characteristics over the total vertical and horizontal flight operating 
regime will require the addition of inlet and nozzle environmental effects 
peculiar to the V/STOL 'system to the conventional engine simulation. 

These effects include the following: , 

• The effects of high angles of attack on inlet recovery and dis- 
tortion during the flight transition between vertical and hori- 
zontal operation. 

• The effects of the above inlet distortion on internal fan and 
compressor performance and stall tolerance, 

• The effects of large exhaust system deflections on nozzle per- 
formance. 

• The effects of flow re-ingestion and associated inlet temperature 
distortion associated with operation in-ground effect or at low 
flight speeds. 

It is expected that the supersonic V/STOL aircraft will operate with 
an open auxiliary inlet during vertical operations and at low forward speeds 
and that the auxiliary inlet characteristics will dominate the level of 
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inlet recovery and distortion. Figure 18 illustrates the effects of 
flight Mach number and mass flow ratio on the inlet recovery character- 
istics of a representative auxiliary inlet. Examination of an inlet 
distortion pattern measured with an open auxiliary inlet indicates that 
the inlet distortion parameter which is defined by: 

ID - (PAV - PMIN)/PAV 

is equal to approximately three times the inlet revovery loss. It has 
been assumed that this relationship would be approximately valid over 
the entire range of flight Mach numbers and inlet mass flow rates such 
that: 

ID ■ 3 (1.0- Inlet Recovery) 

Figure 19 illustrates the effects of this inlet distortion on the fan 
mass flow, efficiency, and stall pressure ratio and Figure 20 contains 
the distortion transfer effects through the front and rear block fans to 
fan discharge. 

Compressor inlet distortion will depend upon the portion of the fan 
exit distortion which is ingested by the compressor and can be approxi- 
mated by the relationship: ^ 

(ID) ■ (ID)^ Compressor Mass Flow 
comp fan — =r — 

^ Fan Mass Flow 

Figure 21 illustrates the resulting inlet distortion effect on compressor 
mass, flow, efficiency, and stall pressure ratio for the baseline VCE/RALS 
engine. 

The engine exhaust nozzles will be subject to significant crosswind 
effects at the low flight speeds associated With the flight transition 
region. The crosswind velocity will be the resultant of the aircraft 
motion and the unsteady environmental air velocity. Figure 22 illustrates 
the effect of the crosswind velocity on the nozzle discharge coefficient. 
Note that the effect is dependent upon the effective nozzle aspect ratio. 
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Corresponding crosswind effects on nozzle thrust coefficient are shown in 
Table 10. Note that these effects are relatively small. 

The inlet and nozzle effects model will be added to the component- 
level engine transient model in order to provide a capability for more 
realistic simulation studies in the hover and low-speed transitional 
region. 


6.2 Aircraft Model 

A simplified dynamic model of a V/STOL aircraft has been modified to 
fit the baseline VCE/RALS propulsion system for use in studying the flight 
transition process. The model uses aircraft stability derivatives obtained 
from Reference 2. Assistance in the development of the dynamic model was 
provided by the McDonnel Douglas corporation. 

The aircraft is controlled aerodynamically during conventional horizontal 
flight by ailerons, rudder and stabilator. Aerodynamic control is augmented 
by thrust control during the hover and low flight speed operating regime. 

The thrust control system assumes a two-engine installation and consists 
of the following: • 

• Fitch control is provided by modulating /he fuel flow to the 
primary and remote augmentors about nominal levels which are 
scheduled throughout the vertical and transition flight regimes. 

The nominal augmentor fuel flow /schedules are designed to produce 
a zero propulsive pitching moment about the aircraft pitch axis. 

• Yaw control is provided by longitudinal vectoring of the ADEN and 
RALS nozzles and by side vectoring of the RALS nozzle. Side 
vectoring of the ADEN nozzle is not available with the current 
ADEN nozzle configuration. 

• Roll control is provided by vertical wing-tip jets supplied by 
continuous compressor bleed air. The flow distribution between 
the two jets is modulated to respond to roll control demands. 

The continuous bleed level is scheduled throughout the vertical 
and transition flight regime to maintain the required level of 
roll control. 

The present model configuration is limited to flight simulations 
initiated at zero forward speed (hover) and terminated at an air speed 
of about 105 knots. Aerodynamic control has not yet been integrated with 
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the model thrust control system and*thrust roll and yaw controls have 
not yet been added. The current model Is, therefore, currently applicable 
only to the simulation of Vertical-to-Horizontal transitions and only to 
the above flight speed limit. Additional limitations on the current model 
include the omission of ram drag and crosswind effects on aircraft motion. 

Preliminary VTH transition trajectories were simulated. The results 
indicated the need for some form of pilot model to control the rate of 
transition from vertical to horizontal flight and the corresponding air- 
craft attitude during the transition. Table 11 contains a simple pilot 
model which was defined and incorporated with the dynamic aircraft model 
for this purpose. Note that this pilot model is applicable only to the 
VTH transition. The pilot model involves the following: 

• Total thrust is held constant at 110% of the aircraft weight 

• Thrust vector angle rotation rate is dependent upon the normal 
acceleration, normal velocity, forward velocity, and pitch attitude 

• Pitch acceleration is dependent upon actual and desired final 
pitch attitude, forward velocity, and the pitch rate 

• Thrust vector angle and pitch rate are obtained by integration 
of the above rates 

The VTH pilot model was utilized in conjunctiot^ with the aircraft model to 
conduct simulation studies of the initial phase of the VTH transition pro- 
cess. The results of these studies are described in Section 7.1. 

6.3 Engine Description 

6.3.1 Engine Components and Gas Flows 

The baseline engine consists of the following turbo-machinery com- 
ponents - a front block fan, a rear block fan, a compressor, a high 
pressure turbine, and a low pressure turbine. Additional components 
include a main burner, two augmentors, and two nozzles. Front block 
fan discharge air can be split between the rear block fan and an outer 
duct which bypasses the air around the second block fan. Second block 
fan discharge air is similarly split between the compressor and an 
inner duct which bypasses the rest of the turbomachinery. The outer 
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bypass from the front block is mixed with the inner bypass air from the 
rear block. The resulting bypass air can be directed to the remote 
burner and nozzle» to the rear mixer for mixing with the hot stream from 
the low pressure turbinei or split between the two systems. Remote 
burner discharge is exhausted through the RALS nozzle and the mixed 
stream from the rear mixer is exhausted through the afterburner and 
ADEN nozzle. Current propulsion control studies have been restricted 
to single bypass operation of the baseline engine with no outer bypass 
flow around the second block fan. 

Gas flow paths are shown schematically in Figure 23. Major components 
are identified by the shaded areas. Stations are identified at the 
boundaries of each component where there is an Interface with another 
component. Thermodynamic accounting is accomplished at each station (or 
node) and Includes accounting for pressure, temperature, enthalpy, gas 
flow, and air flow. In addition to the major gas flows, there are secondary 
flows for fuel flow, cooling air from the compressor to the hot engine 
parts, and customer bleed flow to airframe systems. All flow paths are 
indicated and labeled. The cooling flows are designated as WACLi where 
i is the station number where the flow is assumed to be returned to the 
primary system. The bleed flows are designated as WBi where i is the 
station where the flow is removed from the engine. The fuel flows are 

designated as WFi where i is the last station before the fuel flow is 

introduced into the primary system. Combustion is assumed to be completed 
before the next station. 

A more conventional description of the station designations can be 
obtained from Figure 25. 

6.3.2 Environmental Inputs and Outputs 

The environmental inputs to the engine from the atmosphere and the 
air frame are: 

F2 Inlet Total Pressure 

T2 Inlet Total Temperature 
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PAMB Ambient Static Pressure assumed to be the back 
pressure on the nozzles 

WB28» WB3 Customer bleed air extracted from the compressor 
for use by air frame systems 

The engine outputs are the two thrust vectors. The rear thrust 
vector is obtained from the ADEN nozzle and is defined by a gross thrust 
magnitude and a longitudinal rotation angle from the horizontal direction 
around the pitch axis. The remote thrust vector is obtained from the 
remote (RALS) nozzle and is defined by a gross thrust vector and two 
orthogonal rotation angles - longitudinal rotation and a transverse 
rotation around the roll axis. A corresponding ram drag vector parallel 
to the aircraft velocity vector has been omitted from the initial modeling 
process. 

6.3.3 Engine and Control System Interfaces 

The control system outputs provide additional Inputs to the engine 
model and its evaluation process. The control outputs include the indivi- 
dual fuel flows » variable stator vane positions » and any variable flow 
areas. The current baseline engine requires the following control outputs 

STP22 Second block fan stator position (angle) 

STP49 Low pressure turbine stator vane angle 

A8 ADEN nozzle throat area 

t ■ ■ 

AE80 Remote flow control valve area 

ASS Remote (RALS) nozzle throat area 

AE16 Rear VABI cold-stream discharge area 

WF36 Primary combustor fuel flow 

WF6 Mixed-flow augment or fuel flow 

WF86 Remote augment or fuel flow 

The transition controller obtains the demands for NF36, A8» A88 (at 
high remote flow rates), and STP22 (at low remote flow rates) from the 
multivariable regulator. All other control outputs are obtained directly 
from open-loop schedules. 
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Internal engine variables sensed by the control system include the 
following: 

N2 The front block fan and low pressure turbine rotor speed 

K25 The rear block fan, compressor, and high pressure turbine 

rotor speed 

M93 The gas Mach number at front block fan discharge 

PS3 The compressor discharge static pressure 

T5 The total temperature at low pressure turbine discharge 

T25 The total temperature at compressor inlet 


N25 and T25 are used to calculate the corrected compressor speed which is 
used for open-loop scheduling of the variable compressor vanes which have 
been assumed to track their nominal schedule and, consequently, have been 
omitted from the transient studies. The remaining sensed variables pro- 
vide inputs for the regulator in determining the closed-loop demands. 

External sensed variables and inputs to the control system include 
the following: 

T2 Engine or front-fan inlet total temperature 

PLA Power Lever Angle 

FCD Flight Control Demand or its equivalent (thrust vector 

angle px thrust ratio demands)' ' 

These inputs define the external environment and the current thrust 
management demands on the engine. 

6,3.4 Operating Modes 

The engine and the model are capable of the following modes of 
operation: 

• Horizontal (H) Mode - The H mode is used for conventional horizontal 
flight and has the remote system sh\»t down. All bypass duct flow is 
mixed with turbine discharge flow in the VABI and the resulting flow 
exhausted through the ADEN. The top sketch in Figure 24 illustrates 
the major gas flow paths for the H mode of operation. 

• Vertical (V) Mode - The V mode is used for vertical takeoff and 
landing. All bypass duct flow is diverted to the remote system 
and discharged through the RALS nozzle. The cold side of the VABI 
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is closed and the turbine discharge flow passes through the hot 
side of the VAlil and is discharged through thd ADEN. The bottom 
sketch in Figure 24 Illustrates the major gas flow paths for the 
V mode of operation. 

• Transition (T) Mode - The T mode is used for flight transition 
between the H and V modes. The bypass duct flow is split between 
the remote system/RALS nozzle and the VABI/ADEN. The T mode is 
further subdivided into T1 and T2 modes as follows: 

T1 - the remote system flow is controlled by the RALS nozzle 

T2 - the remote system flow is controlled by the Remote Flow 
Control Valve (RFCV) 

The T1 mode is used for high remote flow rates (the initial phase of 
the V to U transition and the final phase of the H to V transition) where 
the RALS nozzle is capable of adequate flow control, The T2 mode Is used 
for low remote flow rates (the final phase of the V to H transition and 
the initial phase of the H to V transition) where the RALS nozzle is set 
at minimum area and the RFCV must be used to control the flow. 

The engine model reflects each of the above modes of oparation and 
uses different iteration logic for each individual mode. Th^i accounting 
sequence for the remote and bypass systems is also switched for each mode. 

6.4 Model Desc'.?i^?>tlon 

< 

The engine model' represents a non-linear model of the baseline engine 
which has been combined with non-linear models of sensors, actuators, aug- 
mentors, and nozzles developed under the Phase I study. The model repre- 
sents the steady-state and transient performance of the engine. The engine 
dynamics, and corresponding time constants, can be classified as follows: 

• Gas Momentum Conservation 

• Gas Mass and Entropy Conservation 

• Rotor Dynamics 

• Heat Storage 

The current function of the engine model is to develop the engine 
control concept and to demonstrate by simulation that its dynamics are 
suitable for the aircraft control and flight transition processes. The 
gas dynamics effects are about an order of magnitude outside of the 


.003 seconds (time constant) 
.003 to .03 seconds 
.25 to 2.5 seconds 
5 seconds and up 
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desired range and have» therefore, been omitted from the model. Con- 
versely, the heat storage effects have been shown to effect mainly 
acceleration and deceleration rates and some long term thrust transients. 
These dynamics interact very little with the control system and have 
also been omitted. The rotor dynamics have a significant effect on the 
control system and have been included in the model. 

The gas flow and thermodynamic model is a quasirsteady-state model 
with each of the major engine components (fan, compressor, burners, tur- 
bines, mixer, nozzles and ducts) represented as discrete elements of 
the model. Mass flow, pressure, temperature, and enthalpy are accounted 
for directly and the corresponding entropy accounting is handled indirectly. 
The model contains the same inputs and outputs as the engine and the 
order of calculation proceeds roughly from engine inlet to discharge. 

6.4.1 Thermodynamic Accounting 

The turbomachinery components are represented as functions of two 
variables for fixed geometry components and a third variable is included 
for variable geometry components. The input variables for the compression 
components are corrected speed, an arbitrary variable whidh defines the 
departure from the minimum loss line, and variable stator position (where 
applicable). The turbine input variables are" ctarrected speed, an energy 
ratio, and variable stator position (where applicable). Each component 
representation provides the component corrected gas flov?, the pressure 
ratio, and the corresponding temperature ratio. Component efficiency 
is implied by the output pressure and temperature ratios but is not 
calculated. 

Each duct pressure loss is calculated as a function of duct corrected 
flow and the V-BI pressures are obtained from the equations of conservation 
of mass, momentum, and enthalpy. Burner enthalpy rise is a product of 
fuel flow and burner efficiency. Burner efficiency is obtained from the 
burner inlet temperature, pressure and corrected airflow. Nozzle flow 
areas are obtained from choked-flow or unchoked relationships depending 
upon the nozzle pressure ratio. Gas enthalpy is obtained as a function 
of gas temperature and fuel-air ratio and gas entropy calculations are 
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avoided by the direct pressure and temperature ratio calculations. 

The component representations were derived from the steady-state 
component operating characteristics of Figures 8 through 12. Since this 
model is Intended to be a precursor for a subs 'nt real-time model for 
piloted simulation, a minimum computing time appro^iih was used. Conven- 
tional table look-up and Interpolation techniques were replaced by a 
polynomial representation approach. A leaps-and-bounds procedure was 
used to identify the best functional form for each component. The 
Individual coefficients were then determined to minimize the mean 
square error. 

Elemental table-like functions are used for some components. Nozzle 
characteristics, for example, have one function for choked flow and a 
second function for unchoked flow. 

6.4.2 Iteration Techniques 

Internal variables which must be known or assumed in order to complete 
the cycle calculation are designated as iteration variables. An equal 
number of dependent variables must also be available which can be calcu- 
lated by two independent physical processes. The differences between 
the two calculation results represent the iteration errors. The iteration 
errors are functions of the iteration variables and an error partials 
matrix can be calculated from a series' of perturbations of the iteration 
variables. The partials matrix represents a linearization of the functional 
relationships at a specific operating point. The iteration procedure must 
adjust the iteration variables to minimize the error vector. This is 
accomplished by the iteration algorithm which computes an adjustment 
vector from the product of an iteration matrix and the error vector and 
then subtracts the adjustment vector from the previous iteration variable 
vector. 

In a non-real time model, the error partials matrix can be calculated 
for each individual operating point and inverted to obtain the iteration 
matrix. In order to achieve a real-time capability, the dynamic model must 
use a predetermined iteration matrix to permit a single pass per calculated 
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time step. The iteration matrix can be a constant matrix or its elements 
can be scheduled as a function of one or more of the engine variables. 

The iteration matrix will bear some relationship to the inverse of the 
error partials matrix. 

There are two central problems to be solved in order to make this type 
of iteration work effectively. First, the error partials matrix must be 
well conditioned. This must be accomplished by a careful selection of 
iteration variables and an even more careful selection of error variables. 
These variables must be selected or designed for linear independence. A 
measure of linear independence is a condition number defined as the ratio 
of the maximum singular value to the minimum singular value of the error 
partials matrix. If the error partials matrix is A, then the singular 
values are the positive square roots of the eingenvalues of the A'^A matrix. 
For the current engine model which has a 5x5 matrix for the H and V modes, 
and a 6x6 matrix for the T mode, the iteration procedure works well only 
when the condition number is less than 20. This is not a hard limit, 
however, since slow convergence has been achieved with condition numbers 
as high as 150. Extremely rapid convergence has been achieved with condition 
numbers below 15. 

The second central problem involves the scaling of the errors or the 

ft 

iteration variable increments so that the engine error matrix is relatively 

constant for the important coefficients. The variables and errors should 

» 

be ordered so that the most important coefficients are on the diagonal. 
Important coefficients should not vary by more than 4:1 and preferably 

4 

less than 2:1. 

The following procedure was used to develop an iteration matrix from 
the available error partials matrices: 

• The maximum value of the diagonal elements over the set of 
operating points (cycle data used in the regression fitting) 
were selected for the matrix diagonal. 

• The minimum value of the off-diagonal elements were selected 
for the matrix off-diagonal terms. Diagonal terms which 
reversed sign over the operating regime or which were less 
than .1% of the magnitude of the diagonal were set to 0. 
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• The matrix was inverted and then all small off-diagonal terms 
of this inverse were set to 0. The resulting matrix is used 
as the iteration matrix. 

An iteration matrix was developed for each individual model operating 
mode (H» V, T1 and T2) . 

6 . 4.3 Model Structure 

The non-linear cransient engine model has been structured for imple- 
mentation on a hybrid computer with 16 bit fractional arithmetic (EAI P-100) 
All variables have been scaled to operate in the range of -1 to +1. A 
detailed description of the model structure is contained in Appendix C 
showing the functional form of each equation. This description uses 
the following nomenclature: 

P Total Pressure 

Ps Static Pressure 

T Total Temperature 

H Total Enthalpy per Pound of Air 

STP Stator Position 

WA Air Plow ^ 

W Gas Flow 

E Iteration Error , » 

Station designations are illustrated in Figure 25 and are used to 
modify the above variable names. An X following the station designation 
indicates an alternate computation of a variable for the purpose of calcu- 
lating an iteration error. An R following the station designation indicates 
a corrected variable. 

7.0 SIMULATION STUDIES 
7.1 Flight Trajectory Simulation 

The simplified pilot and aircraft models of Section 6.2 were used to 
conduct preliminary simulation studies of the initial phase of the Vertical- 
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to~Horizontal transition process in order to determine thrust require- 
ments during transition. 

Figure 26 contains the results of a VTH transition simulation from 
hover to a final airspeed of 105 knots for a final pitch angle of 2.9®. 

The acceleration required 17 seconds to reach the final velocity and 
resulted in an essentially zero incidence angle for the final 7 seconds 
of the trajectory. RALS and ADEN thrust requirements to fly this tra- 
jectory are illustrated in Figure 27 along with the corresponding thrust 
vector angle. The intent of the thrust vector control was to schedule 
the nominal RALS and ADEN thrust split to maintain a zero propulsive 
pitching moment throughout the transition process. This would be achieved 
by scheduling the flow split between the primary and remote systems and 
corresponding nominal RALS and ADEN augmentor fuel flows as a function of 
Thrust Vector Angle Demand (TVA) . The Pitch Rate (PR) demand would impose 
additional RALS and ADEN fuel flow/thrust modulation requirements which 
should be within the - 12% modulation capability of the current baseline 
engine. As indicated in Figure 27, the actual thrust modulation require?- 
ment (FGMOD/FG) exceeded the - 12% baseline capability near the tail-end 
of the trajectory. 

Examination of the aircraft model and the aerodynamic derivatives 
or the pitch acceleration with respect to the forward velocity revealed 
the source of the problem. Figure 28 illustrates the assumed variation 

J 

of this parameter with forward velocity. The current aircraft model 
aerodynamics are limited to small perturbations about the initial air- 
speed and, consequently, can equate the derivative of the pitch accelera- 
tion (^PITCH'VdV) to ( 3 pi'TCH"/XIVEL) where XIVEL is the difference 
between instantaneous and initial forward velocity. This equality holds 
for small forward velocities but not over the velocity! range of the VTH 
transition process. Further development of the aircraft model will be 
required to resolve this problem so that a more realistic assessment of 
propulsion modulation requirements can be obtained. 
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Figure 28. Aircraft Pitch Acceleration Due To Forward Motion 





7.2 Engine Transient Simulation Results 


The Phase II hybrid computer simulation studies consisted of the 
incorporation of two major modifications to the Phase I simulation 
program: 


• The incorporation of a mode-switching digital transition 
regulator that would handle the V, T, and H operating modes 
along with the VTH and HTV transition paths. 

• The incorporation of an iterative non-linear transient model 
of the baseline engine which represented operation along the 
above transition paths. 


The initial modification to the Phase I simulation involved the 
replacement of the analog engine controller by its digital equivalent. 
The digital controller was then used with the Phase I Linear Engine 
Model and the non-linear actuator, augmentor, and nozzle models to repro 
duce the Phase I simulation results. This modification was successfully 
checked-out in the Takeoff Mode and the digital controller coding was 
used as the basis for coding the transition regulator. 


The next modification involved the replacement of the current digital 
controller by the digital transition controller. The digital controller 
required the following inputs: 


• TRNDIR - This input specified the type of flight transition path. 

A negative value of TRNDIR resulted in the selection of the VTH 
transition and a positive value the selection of the HTV path. 

• PCD - This input represented the Flight Control Demand and defined 
a specific operating point along the specific VTK or HTV path 
defined by TRNDIR. Specific controller proportional and integral 
gains were obtained from the built-in gain schedules corresponding 
to the above inputs. 

• ZEROIC - This third external input permitted the initialization 
of the digital integration process. 

This version of the controller was restricted in that it did not 
permit switching bet’^s^jan single and four-loop controllers along the HTV 
transition and it did not permit reassignment of manipulated variables 
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along either transition path. The digital transition controller was 
then checked-out with the Linear Engine Model for operation at Takeoff 
and Landing. The Takeoff check-out was successful but difficulties were 
encountered at Landing. Program priorities did not permit the effort 
necessary to resolve these problems at this time. 

The Phase II non-linear engine model was then checked-out on the 
hybrid computer as a separate program. Static check-outs (no dynamics) 
were calculated for two VTH operating points (T/0 and MH) and for three 
HTV operating points (FID, RALS, and L) . A number of small problems 
were Identified and corrected. 

The digital controller was then modified to include mode switching 
logic as well as the logic needed for the reassignment of the variables 
ASS and STP22. Note that ASS Is used as a closed-loop manipulated variable 
at high remote flow rates and Is replaced as a manipulated variable by 
STP22 at low flow rates. The coding in this latest version of the con- 
troller contained MIM/MAX statements to limit the outputs of the digital 
integrators as well as numerous look-up tables for establishing upper and 
lower limits on the variables and for providing values for the corresponding 
open- loop schedules. Coding was also added to accommodate a noisy FGD 
signal In order to avoid erroneous manipulated variable or mode switching. 

A special program was developed for exercising the many logic paths 
availably In the controller. 

The revised digital controller was then combined with the non-linear 
engine model In efder to titaln the final Phase II engine and control model. 
An initial attempt was -made to simulate the total VTH transition from 
Takeoff (T/0) to Maximuit Horizontal (MH) . This run was unsuccessful and 
was abandoned. The linear simulation was utilized to generate open-loop 
response data that could be compared with the open-loop hybrid simulation 
data in order to Identify specific problems that could resolve the overall 
simulation problem. The integral portion of the proportional/integral 
controller was disabled in both simulations in order to obtain steady state 
rfs^dnsest The magnitudes and signs of the linear open-loop responses were 
then cbSpared with those &£ the hybrid simulation. This procedure identified 
a nulber of problems and led to modifications of the feed-forward gains and 
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of the nozzle area schedules. (XA8 and XA8B). These changes produced 
hybrid simulation responses that were consistent with the linear 
simulation results. 

The closed-loop hybrid simulation of the VTH transition was tried 
again and was found to work well at Takeoff* Additional problems occurred 
as FGD was varied to move the simulation towards Maximum Horizontal. An 
attempt was made to run at MH but was also unsuccessful. The data from 
these experiments was, however, used to modify the non-linear engine model 
to permit successful simulation runs in the Initial Condition Mode over 
the entire VTH transition path from Takeoff to Maximum Horizontal. This 
simulation could, however, be run in the Operate Mode only at the Takeoff 
point. No attempt was made to check-out the simulation over the HTV 
transition path which was substantially more complex with requirements 
for mode switching as well as reassignment of manipulated variables. 

Figure 29 summarizes the results of the current hybrid simulation 

at Takeoff. The closed-loop control variables were XN2, T5, PS3, and 

XM93 and the corresponding closed-loop manipulated variables were A8, 

WF36, STP49, and A88. The open-loop manipulated variables included STP22 

* 

which received its scheduled value from the controller and AE16 and AE80 
which received their scheduled values directly from table look-up. Note, 
however, that all three open-loop variables were functions of FCD. The 
responses shown are due to a step input of the fan speed (XN2) demand 
and its corresponding interactions with the T5, PS3, and XM93 loops. 

This data from the non-llnear mode switching hybrid simulation compared 
favorably with the prior Phase I results using an elementary multi-variable 
controller and the Linear Engine model. The fan speed dynamics were 
sluggish and asymmetrical. The slow d}mamlcs resulted from Improper time 
scaling of the digital Integration whereas the asymmetry was most likely 
the result of the Increased non-linearity of the simulation. 

8.0 CONDITION MONITORING 

Traditionally, engine condition monitoring systems have evolved in 
parallel with the engine design with little forethought being given to 
objectives or to the ability of the system to satisfy these objectives. 
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Figure 29. Hybrid Simulation Response To Fan Speed Demand at Takeoff 
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The approach has been Co select InstrumenCaClon for clie engine and Chen 
(via analysis and CesC/field experience) Co determine what, if anything, 

Che instrumentation can do to further condition monitoring objectives. 

In spite of this development approach, condition monitoring systems have 
performed useful functions on alrcraft/englne systems. This has been 
possible primarily because engine designs have changed slowly, and lessons 
learned from previous experience have been successfully incorporated into 
new condition monitoring designs. 

The V/STOL aircraft, however, is a radical departure from preceding 
aircraft designs and, consequently, it is expected that a condition moni- 
toring system will serve new roles in V/STOL. To effectively fulfill 
these new roles, the condition monitoring system must be "designed*' rather 
than evolved. 

The specific condition monitoring tasks for the current phase of this 
investigation are to formulate a design strategy and to develop analytical 
tools required to support the design plan. The resulting design strategy 
is described in Section 8.1, and analytical tools are described in Sections 
8.2 - 8.4. It is expected that the design plan will be carried out in 
future phases of the V/STOL development program. 


8.1 Design Plan 

Briefly, General Electric's design strategy for a V/STOL condition 
monitoring system is: 

1. Define objectives for the condition monitoring system, initially 
in terms of capabilities or activities, but ultimately as one or 
more quantitative parameters (denoted objective functions) which 
can be estimated from measured data. Associated with each of 
these objective functions would be an accuracy or repeatability 
needed to satisfy the objective. 

2. Develop a list of possible condition monitoring measurements 
which could be made in support of these objectives. For each 
of the potential measurements, realistic estimates of their 
accuracies and/or repeatabilities would be obtained. 

3. Determine a list of possible engine/control hardware changes 
that could influence the ability of the aircraft to achieve its 
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mission or which could affect the measurements and thereby 
jeopardize the condition monitoring objectives (examples 
are engine component degradation or control schedule mis- 
rigging) . For each of these potential influences, a reason- 
able range for the probable variation would be estimated. 

4. Augment the V/STOL cycle model so that it can compute the 
condition monitoring objective functions and the condition 
monitoring measurements, and so that ffSch of the possible 
engine/control hardware changes is modeled in terms of its 
effect on the cycle. This augmented cycle model can then 

be used to generate linear influence coefficients (derivatives) 
which define the effects of each of the hardware parameters 
on the measurements and the objective functions. 

5. Develop a general algorithm which can be used to estimate the 
objective functions from any arbitrary subset of the pros- 
pective measurements. This algorithm can be used to identify 
the objective function accuracies/repeatabilities based on a 
given instrumentation complement. 

6. Evaluate the ability of selected subsets of measurements to 
achieve the original condition monitoring objectives and, 

on this basis, select one or more instrumentation complements 
which achieve the objectives and are potentially cost effective. 
If specific objectives are not achieveable, then an analysis 
can be performed to determine whether improved measurement 
accuracy can lead to the ability to accomplish the goal, and 
if so, what level of improvement is required. 


After one or more condition monitoring designs have been developed 
using the full cycle deck, a stu4y would be conducted to estimate the 
computer resources needed to implement these systems in service. This 
study would focus on model simplification since this is the area of 
greatest potential payoff. As a result of this analysis the compatibility 
of the condition monitoring system with the control computer would be 
determined. 

The following paragraphs elaborate on these design steps and indicate 
the scope of each effort. 

8.1.1 Identification of Condition Monitoring Objectives 

Typical objectives for an engine condition monitoring system include 
detection of limit exceedances, go/no-go decisions for ensuing flights, 
recognition of sudden shifts in performances, extrapolation of performance 
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trends to plan scheduled maintenance, modular fault isolation to permit 
more cost effective maintenance, identification of out-of-tolerance 
control schedules together with indication of required trim adjustment, 
etc. For the V/STOL aircraft, the go/no-go decision, takes on added 
complexity because of the transition from vertical to horizontal flight. 
This maneuver is not present for standard aircraft and will require careful 
review to choose proper condition monitoring objectives. For example, 
the pilot must be certain that a stall will not be encountered during 
this transition and, consequently, compression stall margins will be 
important objective functions for the condition monitoring system. 

8.1.2 Condition Monitoring Measurements 

A second step in the design process is the selection of a candidate 
list of possible condition monitoring measurements. This candidate list 
will include temperatures and pressures at key locations within the engine, 
rotor speeds, fuel flows, static pressure taps at locations where airflow 
is calculatable, variable geometry position settings, etc. For each of 
these prospective measurements, a model must be developed relating the 
measurement to the cycle representation of the engine so that the cycle 
deck will be capable of predicting the measurement. The accuracy and 
repeatability of the measurement must also be estimated. 

8.1.3 Engine/Control Hardware Changes 

Yet another list to be compiled as part of the condition monitoring 
design process is a list of the engine and control hardware changes that 
can occur to alter the performance of the aircraft and/or the objective 
functions. Somo examples of items to be included on this list are com- 
ponent efficiency degradations; compressor pumping capacity changes; 
mis-rigged control schedules; parasitic flow changes; pressure losses 
due to turning, mixing, etc.; stall line degradation; re-ingestion, etc. 

The effort requires the identification of as many of these hardware pro- 
blems as possible whether they represent long-term changes (such as the 
component efficiencies) or short-term fluctuations (such as re-ingestion). 
These hardware effects must also be incorporated into the engine cycle 
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deck so that their effects on the measurements and the objective functions 
can be determined^ It is also necessary to estimate reasonable limits 
to the hardware variations for subsequent analyses. 

8.1.4 Generalized Algorithm Development 

The algorithm is the tool that will be used to estimate the objective 
functions from the measurements in the presence of hardware variations (and 
measurement errors). Our choice for an algorithm is a maximum likelihood 
estimator such as that described in References 3 and 4. In order to develop 
a specific algorithm, it is necessary to perform the following steps: 

1. Select the measurements that will be used to calculate the 
objective functions (this will be a subset of the possible 
measurements) . 

2. Select the hardware changes (a subset of those present) that , 
are to be estimajied from the available measurements. These 
are analogous to state variables of References 3 and 4. 

3. Use the cycle deck to produce a linear model which relates the 
measurements and tha objective functions to the hardware changes. 

4. Using the weighted least squares approach defined in References 
3 and 4, a maximum likelihood algorithm is thus produced which 
will provide a best estimate of the objective functions given 
the available measurements and the' selected state variables. 

In order to produce the general algorithm that is needed to evaluate 
several proposed instrumentation sets, all that is needed is the generation 
of a complete linear model for all measurements and objective functions. A 
specific algorithm is then simulated by substituting nominal values for 
those state variables that are not to be estimated and using only the 
measurements that are to be available. 

8.1.5 Evaluation of Instrumentation Sets 

In order to evaluate a particular instrumentation set, the procedure 
is to select the set of state variables to be estimated, and then use the 
associated algorithm to perform an accuracy analysis for the objective 
functions. For any specific objective function, there will be some 
particular subset of state variables to be estimated from the measurements 
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that will yield a best accuracy for that objective function. In other 
words, the best accuracy for the objective function will be attained 
when some of the state variables are estimated from the data while others 
are assumed to hold their nominal value. Thus a search scheme will be 
established to choose state variables that will provide a good estimate 
of the objective functions (since there are several objective functions, 
some trade-offs may be desirable) . 

Because of the anticipated scope of this effort, it will not be 
feasible to evaluate all possible instrumentation combinations. Thus 
certain basic sets will be selected for analysis which will Include all 
of the measurements which are expected to be in the basic engine parts 
list, plus some additional measurements that might be condition monitor- 
ing items. These selected sets will be evaluated and in addition the 
impact of adding/deleting certain additional measurements to the sets 
will be evaluated. In the present effort, those objective functions 
which relate to the transition from vertical to horizontal flight will 
be emphasized. 

8.1.6 Other Condition Monitoring Tasks 

Having selected some representative instrumentation sets for further 
consideration, additional analysis will be performed to evaluate the 
computer resources needed to Incorporate the algorithm into an onboard 
system. The algorithm described above expects all of the data to be 
obtained at a single test condition (altitude, Mach number, type day, 
humidity, power setting). In practice, the data will be obtained at 
many different test conditions. Thus, a model of the engine is needed 
as part of the condition monitoring software. The size of the software 
is largely dependent upon the required complexity of this model (to pro- 
vide the needed accuracy). Thus, if more instrumentation error can be 
tolerated (with the objectives still being met), a simpler, less exact 
model is possible. The simplification can pay off in both computer size 
and processor resources. For those functions which logically apply to 
an onboard system (for example, module fault isolation is more logically 
a ground based activity whereas go/no-go is an onboard activity) an 



analysis will be performed to estimate the onboard computer resources 
required and specifically to evaluate compatibility with the control 
computer. 

8 . 2 Maximum Likelihood Estimation 

The primary problem of V/STOL engine condition monitoring is that 
variables such as thrust* compressor stall margin, and attitude control 
ability, as well as hardware performance parameters such as efficiency, 
flow capacity, and leakages cannot be directly measured during flight. 
However, properties of the engine such as pressure, temperature, shaft 
speed and stator position can be measured. Thus values of parameters 
which cannot be measured must be estimated by means of an analytical 
model of the system and available measurements. 

An additional complexity is that measurements are difficult to 
interpret due to engine fluctuations and measurement error or noise. 

To solve this problem a statistical approach is needed. 

Several potential approaches to- solving this problem were considered. 

Of these, maximum likelihood estimation was selected as the most promising, 
primarily because it makes use of the greatest amount of information. In 
addition to considering the measurements themselves, this approach can 
integrate the following factors into the solution; 

1) variances and covariances of the measurement error 

2) nominal, "expected" levels for the hardware performance 

3) variances and covariances of the hardware performance parameters 

Given this information, maximum likelihood estimation distributes the 
original discrepancy (between the measured values and expectations based 
on moninal hardware levels) between measurement errors and deduced deviations 
of the hardware from nominal expectation. The resulting solution is the 
most probable solution given all of the available information. 
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8 . 2 1 1 Measurement and Objective Function Models 


In the most general cases, the maximum likelihood (estimation problem 
is non-linear requiring the application of a numerical optimization pro- 
cedure to determine the solution. However, the linearized version, which 
assumes Gaussian probability distributions and requires the use of a 
linear engine model, is computationally simple and convenient. The linear 
requirement is readily accomplished by centering the analysis at the 
expected performance level, and the solution accuracy is not significantly 
impaired if actual hardware deviations are sufficiently small. 

The linear model for the measured variables, z, may be written as 

z » Hx + V (8.1) 

and a model for the objective functions is 

y » Gx • (8.2) 

where 

z is an mxl measurement vector 
H is an mxn model matrix 
X is an nxl state variable vector 
V is an mxl measurement error vector 
y is a pxl objective function vector 
G is a pxn model matrix 

The vector, x, represents the hardware condition of the engine and Includes 
parameters such as component efficiencies, flow capacities and leakages. The 
vector, y, represents objective function variables such as thrust and com- 
pressor stall margin which cannot be directly measured. The G and H 
matrices are linear approximations to the engine model at the expected 
flight condition. To preserve the linearity z, x and y should, be con- 
sidered as representing deviations from the expected flight conditions. 

The measurement error vector, v, is assumed to be a zero mean Gaussian 
error whose variability is represented by the (mxm) covariance matrix 

R = cov(v) (8.3) 
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Finally, the state vector, x, even though constant during measurement, 
is assumed to belong to a population of known statistical properties. 
Specifically, x is assumed to be Gaussian with mean x and (nxn) covariance 
matrix 

M “ cov(x) (8.4) 


8.2.2 Estimation Model 

Given the above definitions, it can be shown (see References 3 and 4) 
that the most probable estimate, x, is given by the expression 

X ■ Kz (8.5) 


where the gain, K, is given by 


and 


T -1 

K » PK R 


“1 T -1 ”1 

P =» (M +H R H) 


( 8 . 6 ) 


(8.7) 


This algorithm is straightforward and the solution is readily obtained 
if H, M, and R are known. After x is calculated, the estimates of both 
the measured variables and the objective functions can be computed from 


and 


Hx 


Gx 


( 8 . 8 ) 


(8.9) 


If V and z are random variables, then x, z and y are random variables 
with expected values of x, Hx and Gx, respectively. It can be shown that 
the covariance of the associated estimate errors are 


, cov(x-x) = P 

cov(z-Hx) «= HPH^ 

T 

cov(y-Gx) = GPG 

It has also been shown in the literature that the maximum likelihood 
estimator is the optional estimation with the covariances given by 
Equations (8.10-12) being the smallest (minimum error variance) of all 
possible linear estimation techniques. 


( 8 . 10 ) 

( 8 . 11 ) 

( 8 . 12 ) 
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8.2.3 Optimal Filter 
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As successive sets of measurements become available for a system, 
the maximum likelihood estimation algorithm can be applied sequentially 
to successively Improve estimates. This process Is called filtering, 
and the recursion equations for the algorithm after the 1th measurement 
became 



(8,13) 

T -1 

(8.14) 

P^H R ^ 

(P^Ii +hV^h)”^ 

(8.15) 

Hx^ 

(8.16) 


It can be shown*that the error covariance, P^, and the filter gain, K^, 
monotonlcally approach zero in the estimation sequence. 


8.2.4 E xample 

Consider the simple gas turbine used to generate net power, P^, as 
illustrated in Figure 30. The expressions for pressure and temperature 
ratio shown in Figure 30 were derived from a straightforward thermodynamic 
analysis assuming one half of the turbine power is required by the com- 
pressor, the working fluid Is a perfect gas, and a net power level of 

P = 2WC T . 
w p o 

A computer subroutine was constructed for Monte Carlo simulation of 
the simple gas turbine model. The compressor and turbine efficiencies 
were considered to be state parameters with true values of 0.845 and 0.895, 
respectively. The pressure and temperature ratios were considered to be 
measured variables with Gaussian measurement noise at the 1 % level. That 
is: 

z(l) = (Pl/Po) + 0.3RN(1) 
z(2) » (T3/To) + 0,3RN(2) 

V7here RN(1) is a standardized norm'll random variable generated by a random 
number subroutine. 
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A filter subroutine was constructed using Equations (8.13-16) 
The H iiiatriK was determined by taking derivatives of the measured 
variables with respect to the state parameters. That is* 


, )2.5 


Sit 


a(I3/Io) 

Si- 


’t "i 


Millio) 


c ' t 


2 -( 1+2 V ) 


The expected values of compressor and turbine efficiency were 
estimated as 0.85 and 0.90, respectively. After substituting these 
values into the derivative expressions, the following H matrix was 
obtained; 


H 


[ 83.851 0 “I 

-3.076 7 . 843 J 


The initial estimate of the state parameters (efficiencies) were assumed 
reliable at the 1 % level. Thus, 


M 


[ l.E-4 0 1 

0 1 .E- 4 J 


The variability of the measurement noise was also specified at the 1% 
level. Thus, 


^ ^ 9.E-2 0 

0 9.E-4_ 


Simulated measurements from the Monte Carlo subroutine were input to the 
filter subroutine for 60 successive flights (estimations), and the results 
are plotted in Figures 31 through 35. The true values for the simple gas 
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turbine fiystem atr® fl©e;S:e<J square sutabolsi, §ni ;fiB>u3Lat#i "fteisured*’ 
values of pressure and temperature are piottsd as dlaffiSnd f^3^0l,§» Pete 
that simulated measurement ueise causes the ^'raMfUfid" valuet t® 
randomly scattered about the constant |rue values itt figprii 31 SFjd 32‘ 

Estimated values are plotted as triangular symbels.> Pet# the 
estimates converge to the true values In Figures 31 through 34 » and the 
estimated values of state parameters (efficiencies) In Figures 33 and 34 
are within 0«1% of the true values after about fifteen estimations. This 
result agrees with the predicted standard deviation of the estimate errors 
(square roots of the trace elements In the computed P covariance matrix) 
plotted in Figure 35. 


8.3 



A significant deficiency of the filter algorithm described in Section 
8.2.3 is that state parameters are assumed to be constant during the 
measurement sequence. Filter estimates can follow small changes In state 
parameters, however, they will lag large changes. In some condition 
monitoring applications the estimation algorithm must be able to follow 
large changes In the values of state parameters. Two methods to adjust 
the optimal filter process in order to Increase its gain are described 
in the following sections. 


8.3.1 Fictitious Process Noise 

The system model described in Section 8.2.1 can be adjusted to include 
random fluctuations of the state parameters. Such fluctuations are called 
process noise and are modeled as: 

\ " ""l-l + 

where w is the nxl process noise vector assumed to have a zero mean Gaussian 
distribution with Variability represented by the (nxn) covariance matrix 

Q = cov(w) (8.17) 


86 



PRESSURE RATIO 


□ ZTRUE 


ORIGfNAt PAdE IS 
OF POOR QUALITY 

O ZMEASURED A ZESTIflATED 



07/15/81 
10 . 948 


30 

FL r GHT 


PLTOO 1 


Figure 31. Pressure Ratio Variations 
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Figur« 35. Predicted Standard Deviation Of The Estimated Errors 
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The addition of process noise caused the optimal filter algorithm to change 
only slightly. That ls» Equation (8.15) becomes 

?! - (8.18) 

where 

These equations Include the Increased variability caused by process 
noise. The covariance matrix, P^, will increase in magnitude causing an 
increase in the gain matrix, K^. Thus, gain can be artificially increased 
by Including fictitious process noise In the filter algorithm even though 
process noise may not be present in the measured system. 

The simple gas turbine model described in Section 8.2.4 was used to 

investigate this technique. The Monte Carlo simulation subroutine was 

modified to include a 1% step change in the true compressor efficiency 

after thirty measurements, but the true value of turbine efficiency was 

held constant. The filter subroutine was modified to Include process 

noise, and four simulation runs were made with fictitious process noise 

—8 "*7 “6 

at values of 0, 10 ,10 , and 10 times the Identity matrix. 

Results for filter estimates of the state parameters are plotted in 
Figures 36 and 37. It can be seen that the estimate with Q = 0 does not 
follow the step jump well, but filter response improves as Q increases 
until the gain becomes so large that the filter estimate becomes adversely 
affected by measurement nclse. Thus, the fictitious noise technique 
involves n, trade-off between improved filter response to changes in the 
state parameters and increased variability of the filter estimates. 

8.3.2 Moving Window 

A disadvantage of the fictitious process noise technique is that 
extra information must be supplied to the filter algorithm prior to 

estimation. An alternative approach which reduces the significance of 
prior information is the moving window concept . 
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Figure 36. Compr««aor Efficiency Variations With Process Noise 
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Figure 37. Turbine Efficiency Variecions With Proceee Noise 
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Measurement error statistics are determined by using a sample of 
the most recent N measurements to estimate the variance of the data. 
That is, after the ith measurement, 



i 

L 

k-i"N+l 


(Zj),/N 





N is chosen large enough to get a reasonable estimate of the true measure- 
ment covariance and small enough that instrumentation and hardware changes 
are negligible in the interval. 


The filter gain and response to state parameter changes can be 
Increased by using a finite memon' modification to the filter algorithm. 
That is, the recursion equation fvor the estimate covariance will Include 
InfOTmation based upon only the most recent N measurements, and Equation 
(8.15) of the filter algorithm is replaced by 

+ H H - H H) , (8.20) 

Thus, if R. is constant, then P. and K. will be constant after the Nth 
i i i 

measurement and will not asymptomatically approach aero as in the optimal 
filter algorithm. 

The simple gas turbine model described in Section 8.2.4 was used to 
investigate this technique. The Monte Carlo simulation subroutine was 
modified to include a 1% step change in the true compressor efficiency 
after thirty measurements, but the true value of turbine efficiency was 
held constant. The filter subroutine was modified to Include the moving 
window equations, and three simulation runs were made with window sizes 
of N ■ 10, 20 and 30 measurements. 

Results for filter estimates of the state parameters are plotted in 
Figures 38 and 39. It can be seen that the estimate with N*=30 does not 
follow the step jump well, but filter response increases as N decreases 
until the gain becomes so large that the filter estimate becomes adversely 
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Figure 38. Compressor Efficiency Vsrlaclons With Moving Window Concept 
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Figure 39. Turbine Efficiency Verieclone With Moving Window Concept 












affected by measurement noise. Thus, both the moving window and the 
fictitious process noise modifications to the optimal filter algorithm 
involve a trade-off between Improved filter response and Increased 
estimate variability. 


8.4 Filter Design and Evaluation 

Two fundamental design problems arise when filter algorithms are 
applied to practical systems: 


1, The truly optimal filter design must include all error 
sources In the system. However, this would place an 
Impossible burden on the available computer resources 
if all possible hardware parameters In a gas turbine 
propulsion system were Included. Thus, state parameters 
with less significant errors must be deleted and a sub- 
optimal filter must be designed which permits condition* 
monitoring objectives to be met. 

2. The optimal filter design requires that exact values of 
error statistics be known. However, these statistics 
are never known exactly, and approximate values must 

be assumed. 

Because of the need to delete state parameters,, and because informa- 
tion about error statistics is not known precisely, analyses should be 
performed to determine the sensitivity of a filter design to differences 
that exist between the filter and one that fits the optimal mold. This 
process Is called sensitivity analysis. It could be accomplished using 
the Monte Carlo simulation technique demonstrated in Sections 8.2.4 and 
8.3; however, a more convenient accurate and economical method using a 
"truth" model and covariance sensitivity algorithm is described in Chapter 
7 of Reference 4. A computer program was developed to investigate this 
design tool and results are presented in the following Sections. 


8.4.1 Truth Model and Covariance Sensitivity 

An equation for covariance of estimate error for the state vector Is 
derived by substituting the basic measurement and estimation model equetions 
(8.1, 8.13, 8.16) into the covariance definition. Thus, 
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- covj^d-nj 11) 

or - (I-K* H) Pi.a,(I-Ki H)'' + K*RK*T 


( 8 . 21 ) 


Equacion (8,21) is the realistic estimate oP the covariance of state 
estimate error for a filter with arbitrary gain. That is, it is valid 
for both optimal and suboptimal filter algorithms. The asterisk super- 

it 

script is placed on the gain symbol» K^, as a reminder of this fact, 

The truth model sensitivity analysis is a very simple concept as 
illustrated in figure 40. A,ssum6d design values for an M-state filter 
model (PJ, R , U are substituted into the filter covariance algorithm. 

In the absence of process noise, for example. Equations (8,14 and 8.15) 
or (8.14 and 8.20) with asterisk superscripts might be used. The 

it 

resulting sequence of filter gains, K^, together with the N-state truth 
model parameters (P^, R, H)» are substituted into the sensitivity covariance 
algorithm, Equation (8.21), to compute the realistic covariance matrix, 

P^, if the design filter were used to process data from the truth model. 

The filter-indicated covariance matrix, P^, will always be optimistic 
and smaller (less accurate) than the realistic "truth" estimate, P^. Thus, 
a true evaluation of a filter design can be made, and the consequences of 
design changes can be accurately determined. 


8.4.2 Examples 

Two applications of the truth model covariance analysis were carried 
out. The simple gas turbine model described in Sections 8.2.4 were assumed 
together with the optimal filter equations (8.13-8.16). The filter design 
assumed both measurerr.ent noise and initial uncertainty in parameter estimates 
to be at the 1% level. 
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Objectives; 

• Reduced State Analysis 

• Sensitivity Analysis 

• Error Budget Analysts 

• Overall Performance Analysis 


Figure 40* Truth Model Sensitivity Analysis Concept 
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The first analysii consisted of a sensitivity analysis of the filter 
design to measurement noise? The resulting standard deviations of the 
estimated errors are shown in Table 12. Case 1 assumes the actual measure- 
ment error to be at the TA leveli and Case 2 assumes actual measurement 
noise at the 0.5% level. The truth- model results show that the standard 
deviation of estimation errors would be almost doubled for Case 1 and 
almost halved for Case 2. I 

The second analysis consisted of a reduced state analysis of the 
filter design. The resulting standard deviation of the estimate errors 
are shown in Table 13. Case 1 assumes that the filter model has the same 
state parameters as the truth model; and the results are identical. 

However* Case 2 assumes no compressor efficiency error in the filter 
model, and Case 3 assumes no turbine efficiency error in the filter 
model. 

Considering the first column of numbers in Table 13* the filter model 
predictions for the reduced state Gases 2 and 3 are smaller and apparently 
better than the results for Case 1. However, the filter predictions for 
Gases 2 and 3 are unrealistic* and should not be compared! Rather, the 

f 

tc^A^lstic truth model predictions in the last column should be compared 
to each other. The truth model predictions for the reduced state Cases 

2 and 3 are larger and* as expected, poorer' estimation models than the 
model for Case 1. This reflects the fact that estimates for Cases 2 and 

3 have larger uncertainty than Case 1* because a source of error was 
deleted In Cases 2 and 3. 

It is expected that the truth model and covariance sensitivity algorithm 
will be a very useful analytical tool in carrying out the design plan defined 
in Section 8.1 of this report. 

9.0 LONG-RANGE TECHNOLOGY PLAN 

A long-range V/STOL propulsion control technology plan was developed under 
the Phase I program In order to Identify data base and technology requirements 
for the design of Integrated aircraft-engine control systems for supersonic 
V/STOL. Each Individual technology program has been refined to focus on 
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TABLE 12 . SENSITIVITY ANALYSIS OF MEASUREMENT NOISE 



STANDARD DEVIATIONS 
ERRORS 

> 

VARIABLE 

FILTER 

MODEL 

PREDICTIONS 

TRUTH 

MODEL 

PREDICTIONS 




COMPRESSOR EFFICIENCY 
TURBINE EFFICIENCY 
PRESSURE RATIO 
TEMPERATURE RATIO 

(CASE 1 : R ■ 4R*) 

3.34E-3 

3.75E-3 

2.80E-1 

2*80E-2 

6.37E-3 

7.06E-3 

5.34E-1 

5.34E«2 

COMPRESSOR EFFICIENCY 
TURBINE EFFICIENCY 
PRESSURE RATIO 
TEMPERATURE RATIO 

(CASE 2: R « R*/4) 

3.34E-3 

3.75E-3 

2.80E-1 

2.80E-2 

1.96E-3 

2.26E-3 

1.64E-1 

1.64E-2 

SPECIFIED PARAMETERS: 

1 ■ 1 

Po* * ^ 0 * (2,2) - l.E-4 

R^ (1,1) ■ 9.E-2 
R* (2,2) - 9.E-4 



ORIGHNAl PAQE m 
OP POOR 0*«''HTY 


TABLE 13. REDUCED STATE ANALYSIS 



STANDARD DEVIATIONS 

wmwrmmm 


FILTER 

TRUTH 


MODEL 

MODEL 

VARIABLE 

PREDICTIONS 

PREDICTIONS 




(CASE li 

FILTER MODEL ■ TRUTH MODEL) 


COMPRESSOR EFFICIENCY 

3.34E-3 

3.34E-3 

TURBINE EFFICIENCY 

3.75E-3 

3.75E-3 

PRESSURE RATIO 

2.80E-1 

2.80E-1 

TEMPERATURE RATIO 

2,80E-2 

2.80E-2 

(CASE 2; 

NO COMPRESSOR EFFICIENCY IN FILTER) 


COMPRESSOR EFFICIENCY 

0. 

l.OOE-2 

TURBINE EFFICIENCY 

3.57E-3 

4.95E-3 

PRESSURE RATIO 

0, 

8.39E-1 

TEMPERATURE RATIO 

2,80E-2 

2.83E-2 

(CASE 3: 

NO TURBINE EFFICIENCY IN FILTER) 


COMPRESSOR EFFICIENCY 

3.18E-3 

4.19E-3 

TURBINE EFFICIENCY 

0. 

l.OOE-2 

PRESSURE RATIO 

2.S7E-1 

3.51E-1 

TEMPERATURE RATIO 

0.98E-2 

7.08E-2 


SPECIFIED PARAMETERS; 
i ■ 1 

(M) - Po* (2,2) - l.E-4 
R* (I,V - 9.E-2 
R’^ (2,2) - 9.E-4 


vc'rk caski to li© ot'iintid fcowaitU the development of an 
©mini and eonttol iyatem model \d\lnh tould be need In a piloted airotaft'* 
engine ilmulatlon pvogtam» i’he ovetaU ptopaw haa been veatvuotiu'ed to 
delete the oontvol cowponenct and infeeitation pvoRtawi whioh wete ppimatlly 
hatdwati deaign ovtented and to add a pfiogvam fov V/StOh engine oondition 
fflonitoting'i fhe following iectione deaovibe the individual pvo|tawe and 
an overall plan whlel\ haa been baaed on an eaaentially conatant level'^of* 
effovt fov each yeaw 


0 % 1 Qontral Reauivementa 

'fhe Qonttol veciulifementi pvoivaw ia oonoetned with tbe eatabliahment 
of engine operating aeUedulea and protection icheduleii w'^th parameter 
eenaing rec^uirement# for implementint tl\e acbedulea> and with corresponding 
dealgn margin reqwirementa for the ielected achedviling concepta% individual 
work elementa include the followingi 


• Yvanai tion Qonti’^ola » The tranaition controls repreaenfe the engine 
protection acUedulee which limit regulator action on large and/or 
rapid throttle tranaientst The tranaition control iv\cl.udea the 
acoel fuel ichedule which limits maximum fuel flow during throttle 
hurata to protect against oompreaaor etall and the decel fuel 
achedule which limita minimum fuel flow during throttle chops to 
prote^f againat comhuator blow-out » It also includea protection 
against over-apeed» temperature,, and preaaure during vSteady-atate 
and tranaient operation% The tranaition control invoXvea parameter 
aenaing recivvirementa and proceaaing reqviirementa for engine fan 
and con 8peed» compreaaor inlet temperaturei compreaaor discharge 
preaaurei and at least one turbine temperaturei 

• Augmentor Mght-Qff Awgmentor light-off control a are required for 
both’”' tl\¥ 'pima^^ and remote augmentorii Kach control inclndea fuel 
metering to individual pilot burner and combustion aonea, flame 
xTetection, and nosale pre-opening or internal ge.ometr>’ modulation 
to minimise internal preaaure apikea. 

« Pov^er Kg>aaement - Power management atudiea must be conducted to 
identiiy the apecific aenaed or computed engine parametera which 
should be acheduled againat power lever angle and thrust vector 
angle detwanda for aetting engine thruat magnitvide vand directioni 
Additioxml atudiea twjtt ideutify the correaponding parameters 
C thruat aettiug parameter, corrected fan apeed, corrected core 
speed, etc 4 for acheduling each individual open and cloaed-loop 
demand achePleii 


lOA 



• PvtXX’-Ranftt Cattteyol Sehadule^ « Sueady»8tat»fc operating icliaduXea 

«iuafe ba dttvalopad ^ w nmuageinenu pai'amtai-St eloaad“* 

Xoop aatvaecl vai*iable8| aud ^pen-loop wauipu.lttt«d vatlabXas acroaa 
cha Cull opai*atii\| anvalopa (aXtltuda, Had\ i\un\b«p» and ait\biii\fc 
Kainp®i.*atui.‘a i*ang®a) and fcht full nparatiing I'anga (powai' Xavae 
angla and tjUvnat vacfeot angXt), Any apaclaX ra«i«4rainanfei fioi* bltad 
nnwpanaabion shonXd also be defilned fov fcba Ulgh bleed opei‘autng 
regime (vevcinal and etanattlonal flight regimes). 

• Stall froceetlon - Any speeial eeqnlrementa for fan and compresaoi* 
a tall pretention during high preaam'e diatortion or temperature 
re«ingeation ahould be identified and potential oontrol aolutiona 
examined % 

« Design Marxina - I'be selected oloaed^loop control configuration and 
power manaiement concept must be combined with avaiiabie component 
quality and deterioration models and control tolerance models to 
establish corresponding design margin requirements on rotor apeedsi 
temperatures* pressures , stall margins* etc. 


9 ’ * Modeling and Simulation 

4he modeling and simulation pfogram will be concerned with refining 
the current fhaee XI engine and aircraft models and with additional simulation 
studies to evaluate iarga throttle transients* special environmental effects, 
and stability characteristics over the full operatiug range. Individual work 
elements lx\clude the following; 


• Knaina Model - I'he current engine model must be expanded to include 
the inlet and noKtle effects model, transition controls, and the 
full-range control schedules, Potential uechnitiues for improving 
engine model computer ruuning time must be identified, bested, and 
evaluated, 


• Aircraft Model - The aircraft aiodei must be refined to include force 
and moment coefficients out to the *3 Mach number regime, a ram drag 
vector, and a capability for initiating transitions at ,3 Mach in 
horisonbai flight idle descent, The existing pilot model must also 
be refined to provide input demands to both the propuliior system 
and the flight control system dviring the translbioual flight regime. 
Complete V'i'H and Hl'V transitions should be simulated in order to 
evaluate the current ?/S'X'0b control concept, 

• Invironmanbal Effects ^ Simulation studies must be conducted to 
determine the effects of inlet pressure distortion and temperatnre 
re-ingestion to the operating characteristics of the engine and 
control system, Requirements for special engine protection control 
logic must be identified. 


lOS 


• System Evaluation An evaluation test plan must be developed 
and carried out to evaluate system response characteristics over 
the major flight operating regimes. These studies must Include 
small and large throttle bursts and chops ^ bodes* and mode 
transitions to Identify any potential operating limitations. 

e Real Time Model - All significant run-time refinements must be 
Integrated Into the expanded engine and control model to produce 
a final system model which would be suitable for piloted simulation 
at NASA-Mes. The final model would be coded* checked-out* and 
transmitted to NASA. 

9.3 Control 

The control logic program will be concerned with the integration of the 
Individual sub-systems Into an overall control system logic design* with the 
Identification of Input-output Interface logic with the aircraft/flight control 
system* and with the Identification of preliminary fault management require- 
ments. Individual work elements Include the following: 


• Control System Model - Control system requirements for transition 
control* full-range operating schedules* and regulator gain schedules 
must be implemented In control logic and Integrated with the regulator 
and feed-forward system logic into an overall control system model. 

Any special engine protection requirements or operating refinements 
^ould be added as required. | 

• STOL Operation - The control system model will be evaluated with 
respect to Short Takeoff and Landing (STOL) operations and any 
special control logic requirements will be Identified* developed* 
and integrated Into the control system model. 

• Interface Logic - Interface logic must be developed for processing 
Input demand and environmental signals from the aircraft/flight 
control system and for transmitting output signals on engine 
operating condition* available control margins* and essential 
flight safety information. 

• Fault ManaRment - Fault management studies will be conducted to 
establish preliminary sensor and actuator FiCA logic and to 
evaluate system effectiveness In the vertical* transition* and 
low flight speed operating region. 

9.4 Condition Monitoring 

A condition monitoring program plan was developed and has been described 
In Section 8.1. The Initial phase of this plan was used as the basis for the 
condition monitoring studies summarized in Section 8.0. Subsequent condition 
monitoring studies are expected to be pursued under a separate program and* 
consequently, have been omitted from the overall technology plan. 
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9.5 Propulsion Control TechnoloRy Plan 


Th« Individual technology programs have been integrated into the 
overall propulsion control technology plan summarized in Figure 41. It 
represents a time-phased program which will lead to a real-time V/STOL 
propulsion and control system model suitable for piloted simulation studies 
by mld-1985. The indicated program maintains the current technical level- 
of-ef£ort through each of the indicated program phases. 

10.0 DISCUSSION OF RESULTS 

Typical propulsion requirements were examined for a maximum gross weight 
transition from vertical takeoff to horizontal accelerated-climb and for a 
minimum gross weight transition from flight idle descent to vertical landing. 

Key operating points were identified along each transition process and were 
used to establish typical steady-state operating requirements. The overall 
thrust vector direction was rotated from the vertical to the horizontal 
direction by longitudinal deflection of the ADEN nozzle, by varying the rear 
VABI area to modulate the remote-to-primary flow ratio, and by modulating 
the RALS augmentor fuel flow to maintain a zero propulsive pitch moment 
about the aircraft. Nominal ADEN augmentor temperature was maintained to 
preserve height and pitch control capabilities. Operation at low remote 
flow rates required the use of a remote flow control valve to hold the RALS 
augmentor inlet Mach number in the stable combustion regime. 

Steady-state and transient partial derivatives were calculated at each 
individual VTH and HTV operating point and were used to conduct control mode 
studies and to develop linear state-space models for the regulator design 
process. Tha mode studies indicated that a four-loop regulator design could, 
be retained throughout most of the flight transition regime. Mode transitions 
were: required from closed-loop A88 control to STP22 control at low remote flow 
rates and to a single-loop regulator for flight idle operation in the horizontal 
flight mode. The K/Q matrix technique (described in the Phase I report) was 
used to develop proportional and integral gain constants for the four-loop 
controller for each individual operating point. Linear regression was used 
to establish overall gain schedules for the VTH and HTV controllers as a 
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Figure 41. Long-Range Prr>pulsion Control Technf^logy Plan. 













fuuution of n prallniiimry confvol demand par^metei*. Simulation studies 

indicated acceptable response cluvracteciatics with relatively little cross- 
coupling effects for the application of simultaneous unit step demands to all 
four loops. 

Mathematical models were developed for representing inlet and nottle 
environmental effectsi aircraft flight dynamics, and engine transient char- 
acteristics. I'he inlet and nostle model provides corrections to nominal fan 
and compressor performance for inlet distortion and to the nostle discharge 
coefficient for deflected nosssle operation. It will be ineorporated into 
the engine model at a later date, live aircraft model was used to simulate 
the initial phase of the VTU transition in order to evaluate the sero pro- 
pulsive pitch concept, Additional refinements to this model are needed to 
add the effects of the ram drag vector, combined aero and propulsive flight 
control, and a capability out to .3 Mach number, the engine modeling 
involved the development of a non-linear component-level regression modal 
of available steady-^state operating characteristics, the non-linear com*^ 
ponent model vv’as combined v>?ith the Phase I non-linear models of the actuators, 
augmentors, and nosales i with a modified version of the Phase X feed-fon<?ard 
system*, and wltVi the digital transition controller and its corresponding gain 
schedules. Subsequent non-linear simulation studios on the hybrid computer 
v^^ere successful in running individual operating points along the Vertical-to- 
Uoriaontal (VTH) transition trajectory in an initiallalng mode but were not 
successful in running transients along the VTH path. Corresponding simulations 
of the Uoriaontal-tQ-VerUical (KTV) transitions were not attempted because of 
the greater complexity of the HTV transition. 

The simulation study results did, however, provide a partial validation 
of many of the design concepts and procedures employed in the development of 
the non-linear engine model and of the transition controller. Successful 
design concepts and procedures included the following; 

• The use of multiple regression techniques for the development of 
an autcmatically scaled non-linear component-level engine model. 

• The use of linear simulation studies for designing a multi-variable 
regulator with sufficient robustness for surviving tlve transition 
from the linear world to the non-linear world. 
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• The use of the linear simulation for generating test responses 
for checking the non-linear simulation. 

• The concepts of mode-switching and reassignment of manipulated 
variables for handling regulator configuration changes. 

The lack of success in simulating any significant transients along 
the VTH (or HTV) transition paths indicates a need for further development 
along the following lines: 

• Better definition and/or refinement of the closed-loop and open-loop 
control schedules along the VTH and HTV transition paths. 

• Better definition and refinement of the feed-forward schedules 
over the VTH and HTV transition flight paths. 

• Further development and orchestration of checkout procedures 
utilizing linear and non-linear time-sharing studies in con- 
junction with the non-linear hybrid simulation studies. 

Preliminary V/$TOL engine condition monitoring studies were initiated 
to establish system design objectives, measurement requirements and evaluation 
criteria, and design techniques and procedures for developing a viable con- 
dition monitoring system. A simple gas turbine cycle was established as an 
analytical test vehicle and was used to evaluate a number of potential 
approaches for maximum likelihood estimation of engine health parameters 
which are not directly measureable (such as thrust and stall margin). 
Comparative results were obtained for a moving window approach and for an 
approach using fictitious process noise which is provided externally. A 
technique for evaluating reduced state models (which could be accommodated 
in the propulsion control system) against a truth model was also examined. 

The long-range technology plan was re-structured and refined to focus 
specifically on propulsion control concepts and capabilities required for 
the development of a high fidelity real-time simulation capability which 
could be used in a piloted aircraft-engine simulation program. The revised 
plan describes individual sub-programs in the areas of control requirements, 
modeling and simulation, and control logic. The overall program assumes that 
the current level-of-effort would be maintained constant, and therefore, 
extends for four more years through the 1984-85 time period. 
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11,0 APPENDICES 


11,1 Appendix A - Individual Regulator Response To Unit Step Demands 

Appendix A contalne time domain plots of the effects of simultaneous 
unit step demands on all inputs for the individual regulator designs of 
Tables 6 and 7. Figures Al<-A6 contain the results for the individual VTH 
regulators » and Figures A7-A11 for the multi-loop HIV regulators. Note 
that) in all cases» multi-loop interactions appear to occur in about the 
first ,25 seconds and that steady-state response is achieved by 2 seconds. 

The following nomenclature has been used in Figures A1 through A22: 


U1 

Fan Speed (PCN2) Demand 


U2 

Turbine Discharge Temperature 

(T5) Demand 

U3 

Compressor Discharge Pressure 

<PS3) Demand 

U4 

Duct Mach Number (XM93) Demand 


Y1 

Fan Speed Response 


Y2 

Turbine Discharge Temperature Response 

Y3 

Compressor Discharge Pressure Response 

n 

Duct Mach Number Response 

* 


11 , 2 Appendix B - Scheduled Regulator Response To Unit Step Demands 

Appendix B contains time domain plots of the effects of similar unit 
step demands for the VTH and HTV regulators based on the linear regression 
fits indicated in Tables 8 and 9. Figures A12-A17 contain the results at 
the individual operating conditions for the VTH transition, and Figures A18- 
A22 for the corresponding HTV operating conditions. These results Indicate 
somewhat greater multi-loop Interactions and, in some cases, longer time 
requirements to achieve steady-state conditions. 
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FIGURE A-1 RESPOfISE TO UNIT STEP DEMANDS AT TAKEOFF 




FIGURE A-2 RESFOfiSE TO UNIT STEP DEMANDS AT. START TRANSITION 
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FIGURE A-3 RESPONSE TO UNIT STEP DEMANDS AT ACTIVATE RFCV 
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FIGURE A-4 RESPOnSE TO UNIT STEP DEMAfiDS AT MID TRANSITION 
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FIGURE .-6 RESPOMSE TO UNIT STEP DEMUNOS AT MAX,E„„ W,rizo,;tal 
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FIGURE A-7 RESPONSE TO UNIT STEP DEMANDS AT 30? RALS FLOW 
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FIGURE A-9 RESPONSE TO UNIT STEP DEMANDS AT MINIMJM RALS (HTV) 
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FIGURE A-10 RESPONSE TO UNIT STEk DEMANDS AT 80% RALS FLOW 
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FIGURE A-11 RESPONSE TO UNIT STEP DEMANDS AT LANDING 



SCHEDULED REGULATOR RESPONSE TO UNIT STEP DEMANDS AT TAKEOFF 





FIGURE A-13 SCHEDULED REGULATOR RESPONSE TO UNIT STEP DEMANDS AT START TRANSITION 
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FIGURE A-16 SCHEDULED REGULATOR RESPONSE TO UNIT STEP DEMANDS AT MINIMUM RALS (VTH) 
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SCHEOULEO REGULATOR RESPONSE TO UNIT STEP DEMANDS AT MAXIMUM HORIZONTAL 
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FIGURE A-18 SCHEDULED REGULATOR RESPONSE TO UNIT STEP OEHAWS AT 30J RALS FLOW 


Y1 = N2 (Output) 
Y2 = T5 (Output) 
Y3 = PS3 (Output) 
Y4 = M93 (Output) 
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FIGIIRF A-19 SCHEFXILED REGULATOR RESPONSE TO UNIT STEP DEMANDS AT 60? RALS FLOW 
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FIGURE A-20 SCHEDULED REGULATOR RESPONSE TO UNIT STEP DEMANDS AT MINIMUM RALS (HTV) 
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FIGURE A-22 SCHEDULED REGULATOR RESPONSE TO UNIT STEP ^DEMANDS AT LANDING 


11.3 Appendix C - Non-Linear Transient Model 
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Front Fan 


Loss Parameter 
Airflow 

Discharge Pressure 
Discharge Temperature 

Second-Block Fan 

Loss Parameter 
Corrected Airflow 
Inlet Pressure 

Discharge Pressure(to core) 

Discharge Temperature 

Discharge Pressure (to 
Bypass Duct) 

Discharge Temperature 

Airflow. 

Front Fan Discharge 
Mach No . 

Discharge Enthalpy 

Compressor 

Loss Parameter 
Inlet Pressure 
Air Flow 

Bypass Duct Airflow 

Discharge Pressure 

Discharge Temperature 

Inlet and Discharge 
Enthalpies 

Discharge Flow 


6H2 is an iteration variable 
W21 • P2//f2*f(N2//f!, GH2) 
P21 • P2^f(N2//r2, GH2) 

T22 » T2*f(N2/4T2, GH2) 


GH22 is an iteration variable 
W22R * f(N25/jT22, GH22, STP22) 

P22 » P21*f(W22R) 

P22X « TT52*W2/W22R 

P23 « P22’^f(N25/4T22, 6H22, STP22) 

T25 - T22*f{N25/rf^. GH22, STP22) 

PI 3 * P23’^(N25//f22, GH22, STP22) 

T14B « T25*f(N25//f22, GH22, STP22) 

W22 * P22*W22R/Jf?2 

M93 * f(W22R) 

H14B « f(T14B) 

GH25 is an iteration variable 
P25 ^ P23/f(W25R) 

U!2& = W25R*P25//T25 

W14B » W22-W25 

P3 - P25^f(N25//r25, GK25) 

T3 = T25*f(N25/jT25, GH25) 

H25 = f(T25) 

H3 = f(T3) 

W3 = k*W25-WB28 

(k reflects cooling air extraction and 
WB28 customer bleed extraction) 
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Compression Power 

Front-Fan 
Second-Block Fsn 

Compressor 

High Pressure Shaft 

Cooling Flow Accountability 

Combustor 
HP Turbine 

LP Turbine 

Frame 


ORlQiNAL PAGE IS 


PW2 • kl*T22*W2-k2*T2*W2 
PW22 * -k3*T22*W22+k4^H25*W25 
+k4*H14B*W14B 

PK25 * k5*H3*W25-k6*H25*W25 
-k7*H3n(B28 
PWCMP - PW22+PW25 


WA36 = W3-k8*W25 
WA41 = WA36+k9^W25 
WA42 * WA41+klO*W25 
WA49 » WA42+kll*W25 
WAS « WA49+kl2*W25 
WAS6 * WA5+kl3*W14B 


Combustor 

Compressor Discharge 
Pressure. 

Fuel -Air Ratio 

Ef f 1 ci ency 

Discharge Enthalpy 

Discharge Temperature 

Discharge Pressure 


PS3 * P3*f(W3‘fr3/P3) 

FAR36 « kl4*WF36/WA36 
E4D36 = f(P3, W3JT3/P3, FAR36) 
H4 * H3+kl5’^FAR36*E4D36 
T4 * f(H4. FAR36) 

P4 • P3’*f(W3jT3/P3. FAR36) 


High Pressure Turbine 
Gas Flow 

Nozzle Enthalpy 
Fuel -Air Ratio 
Unbalanced Power Factor 
Power 

Rotor Acceleration 
Energy Ratio 
Flow Function 
Alternate Inlet Pressure 
Exit Pressure 


W4 = WA36+kl6*WF36 

W41 » W4+kl7*W25 

H41 = (WA36*H4+kl8*W25*H3)/WA41 

FAR41 = kl9*WF36/WA41 

PWX4Q is an iteration variable 

PW4 * PWCMP (1+PWX4Q) 

DN25 = PWX4Q*PWCMP/N25 
H41R = k20*PW4/(T41*W41) 

W4R = f(N25/ T41, H41R) 

P4X « T4*W4/W4R) 

P42 * P4X*f(H41R, N25/ R41) 


ORIGINAL PAGE IS 
OF POOR QUALITY 


Exit Gas Flow 
Exit Enthalpy 
Exit Temperature 

Low Pressure Turbine 

Gas Flow 

Cooling Flow Enthalpy 
Gas Enthalpy 
Temperature 

Unbalanced Power Factor 
Rotor Acceleration 
Discharge Enthalpy 

Energy Ratio 
Flow Function 
Inlet Pressure 
Discharge Airflow 
Discharge Temperature 

VABI 

LPT Discharge Flow 

Hot Stream Discharge 
Pressure 

Gas Flow 

Discharge Enthalpy 

Discharge Temperature 

Duct Discharge Pressure 

Cold Stream Area 

Hot Stream Area 

Hot Stream Flow 

Hot Stream Static 
Pressure 


W42 * K41+k21*W25 

H42 - (WA41*H41-PW4+k22*W25*H3)/WA42 
T42 - f(H42, WF36/WA42) 

W49 - W41+k23*W25 

HCL49 = k24*H3+k25*H25 

H49 * (WA42*H42+k26-W25*HCC49)/WA49 

T49 = f(H49, WF36/WA49) 

PWX48Q is an iteration variable 
DN2 * (PWX48Q*PW2)/N2 
H5 * (WA49*H49-PW48+k27*W25* 
HCL49)/WA5 

H49R * PW48/{T49*W49) 

W42R = f(N2//T49, H49R, STP49) 

P42X ' 42*P42/W42R 
W5 « W49+k28*W25 
T5 ■ f(H5, WF36/WA5) 

W5R = W5*JT5/P5 

P56 = P5*f(W5R) 

WE6 = W5+k29*W14B 

H56 = (WA5*H5+k30*W14B*H14B)/WA56 

T56 = f(H56. WF36/WA56) 

PI 6 = P13*f(W13R) 

AE16 is an iteration variable 

AE56 = k31-AE16 

W56RQA « W56* T56/(AE56*P56) 

PS56 * P56-f(W56RQA) 
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Gas Flow 

Remote Flow Ratio 
Remote Airflow 
Cold Stream Flow 
Mixed Air Flow 
Mixed Enthalpy 
Nozzle Gas Flow 
Fuel -Air Ratio 
Mixed Gas Flow 
Nozzle Entrance Flow 
Cold Stream Flow 


W14Z » k32^41B-k33*WA56 

WBL16Q is an iteration variable 

W81 * W142*WBL16Q 

Wli6 s W14Z-W81 

WA58 * WA56+WA16 

H58 * (WA56*H56+W16*H14B)/WA58 

WA7 » WA58+k34*WA56 

FAR7 » k34*(WF6+WF36)/WA7 

WS8 « WA58+K35*WF36 

W7 * W5B+k36*WA56+k37*WF6 

W16RQA = W16*ffT4/(AE16*P16) 


Cold Stream Static 
Pressure 

Mixing Plane Momentum 

Mixed Velocity 
Mixed Temperature 
Mixed Static Temperature 
Mixed Pressure 

Afterburner 

Fuel -Air Ratio Added 
Efficiency 
Discharge Enthalpy 

ADEN Nozzle 

Entrance Enthalpy 

Temperature 

Pressure 

Gross Thrust 


PS16 = P16*f(W16RQA) 

WV16 W16*f(T14B, W16, W16RqA)/(P16*AE16 
WV58 * W56*f(T56, W56, W56RQA)/(P56*AE56) 
V58 » (k38*WV56+WV16)/W58 
T58 • f(H58. WF36/WA58) 

TS58 - T48-k39*V58 
P58 « PS56*f(V58, TS58) 


FAR68 = WF6/WA58 

E68D6 - f(P58. FAR68, T58) 

H68 * H58+k40*E68D6*FAR68 


H7 = (WA58*H68+WA56*H14B)/WA7 
T7 = f{H7, FAR7) 

P7 = P58*f(V58, FAR68) 

P7X = f(W7*w'TT/A8, PAMB) 

F69 = W7*f(P7, T7) 


RALS (Calculated only when W81 is greater than zero) 


Alternate Remote Flow 
RFCV Flow 

RFCV Discharge Pressure 


W81X k41*AE80*PT6/JT14B 
W80RQA * W81*iTT^/{P16*AE80) 
P86 * P16^f(W80RQA. AE80) 




Burner Fuel -Air Ratio 
Burner Exit Pressure 
Burner Exit Enthalpy 
Burner Exit Temperature 
Nozzle Gas Flow 
Alternate Nozzle Pressure 
Gross Thrust 

Iteration Errors 

E4 * P4-P4X 
E42 - P42-P42X 
E22 * P22-P22X 
E7 « P7-P7X 
E87 » P87-P87X 
El 6 ■ PS16-PS56 

j 

E81 * W81-W81X 

ESI 87 - MIN(E87, -E81) 


FAR86 ■ k42*WF86/W81 

P87 « P86*f(W81*JrT4B/P86. FAR86) 

H87 * H14B+k43*FAR86 

T87 * f(H87, FAR86) 

W87 - W8Uk44*WF86 

P87X • f(W87*^/A88. PAMB) 

FG89 * W87*f(P87. T87) 


V, H, and T Modes 
V» H, and T Modes 
V, H. and T Modes 
V, H, and T Modes 
Used In E8187 
H and T Modes 

Used In E8187 
V and T Modes 
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11, A Appendix D - Nomenclature 


Symbol 

Description 

Engr, Units 

ADEN 

Augmented Deflected Exhaust Nozzle 


ADL 

ADEN Deflection Angle (longitudinal) 

degrees 

AR 

Nozzle Aspect Ratio 


A/RFCV 

Activate Remote Flow Control Valve (VTH Transition) 


A8 

ADEN Throat Area 

inches'^ 

A27 

VABI Cold Stream Area 

inches'^ 

ASS 

RALS Nozzle Throat Area 

inches'^ 

all-a44 

Regulator Proportional Gain Coefficients 


bll-b44 

Regulator Integral Gain Coefficients 


BNG 

RALS Nozzle Deflection Angle (Transverse) 

degrees 

Cd 

Nozzle Discharge Coefficient 


CFGR 

Nozzle Thrust Coefficient 


Cp 

Specific Heat at Constant Pressure 

BTU/lb/"R 

ETA . 

Fan Efficiency 


FCD 

Flight Control Demand Parameter 


FG 

Gross Thrust 

lbs. 

FG8 

ADEN Gross Thrust 

lbs. 

FGse 

RALS Nozzle Gross Thrust 

lbs. 

F689 

RALS Nozzle Gross Thrust . 

lbs. 

FGT 

Total Gross Thrust 

lbs. 

FID 

Flight Idle Descent (HTV Transition) 


FN 

Net Thrust 

lbs. 

PRAM 

Ram Drag 

lbs. 

G 

Model Matrix (Objective Function) 


H 

Model Matrix (Measured Variables) 


HTV 

Horizontal-to-Vertical Transition 


H41R 

Low Pressure Turbine Energy Ratio (AH/T) 

BTU/lb/“R 

H49R 

High Pressure Turbine Energy Ratio (ah/T) 

BTU/lb/®R 

ID 

Inlet Distortion Level 


K 

Gain 


'Q Matrix 

Regulator Design Procedure 



Symbol 

L 

M 

M 

MH 

M/R 

M/T 

N 

P 

PAV 

PCN2 

PCN25 

PIC 

PITCH 

PLA 

PMIN 

PR 

PREl • 

PRE2 

PS3 

Pt 

PTj/Po 

P/P 

Po 

Pi 

Pw 

P2 

P3Q25 

P13Q25 

P93Q2 

Q 

R 

RALS 

RALS 

RDL 
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Description 

Landing (HTV Transition) 

Flight Mach Number 

State Vector Covariance Matrix 

Maximum Horizontal (VTH Transition) 

Minimum RALS Operation (VTH Transition) 
Mid-Transition (VTH Transition) 

Number of Measurements 
Filter Matrix 

Average Inlet Total Pressure 
Per Cent Fan Speed 
Per Cent Core Speed 
Proportional /Integral Controller 
Pitch Attitude Angle 
Power Lever Angle 
Minimum Inlet Total Pressure 
Pitch Rate Demand 

First Precompensator Matrix >' 

Second Precompensator Matrix 
Compressor Discharge Static Pressure 
Total Pressure 

Nozzle Total -tO“Static Pressure Ratio 

Fan or Compressor Pressure Ratio 

Ambient Pressure 

Compressor Discharge Pressure 

Net Power Output 

Inlet Total Pressure 

Compressor Pressure Ratio 

Rear-Block Fan Pressure Ratio 

Front-Block Fan Pressure Ratio 

Process Noise Covariance Matrix 

Measurement Error Covariance Matrix 

Remote Augmented Lift System 

Minimum RALS Augmentation (HTV Transition) 

RALS Deflection Angle (Longitudinal) 


Engr, Units 


psia 

% 


degrees 

degrees 

psia 

deg/ sec. 


psia 

psia 


psia 

psia 

hp 

psia 


degrees 
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Symbol 

Description 

Engr. Units 

RDT 

RALS Deflection Angle (Transverse) 

degrees 

RFCV 

Rr(Ote Flow Control Valve 


SFC 

Specific Fuel Consumption 

Ibs/hr/lb 

SM2 

Front-Block Fan Stall Margin 

% 

SM22 

Rear-Block Fan Stall Margin 

A' 

K 

SM25 

Compressor Stall Margin 

% 

STP22 

Rear-Block Fan Stator Position 

degrees 

STP49 

Low Pressure Turbine Stator Position 

degrees 

STOL 

Short Takeoff and Landing 


SYS 

Linear Engine Model 


S/T 

Start Transition (VTH Transition) 


t 

Time 

seconds 

TKOA 

Linear Engine Model 


TKOIF 

Sensor Matrix 


TKOIP 

Actuator Matrix 


TK02P 

Controller Matrix 


TR 

Thrust Ratio Demand 

■ 

TVA 

Thrust Vector Angle Demand 


T/0 

Vertical Takeoff (VTH Transition) 


AT 

Augmentor Temperature Rise 

®R 

To 

Ambient Temperature 

“R 

T2 

Fan Inlet Total Temperature 

®R 

T3 

Compressor Discharge Temperature 

"R 

T41 

High Pressure Turbine Inlet Temperature 

“R 

T48 

High Pressure Turbine Discharge Temperature 

“R 

T49 

Low Pressure Turbine Inlet Temperature 

®R 

T5 

Low Pressure Turbine Discharge Temperature 

®R 

V 

Flight Velocity 

meters/sec. 

V 

Measurement Error Vector 


VABI 

Variable Area Bypass Injector 


VCE 

Variable Cycle Engine 


V/STOL 

Vertical and Short Takeoff and Landing 


VTH 

Vertical-to-Horizontal Transition 


WF36 

Primary Fuel Flow 

Ibs/hr, 

WF6 

Mixed-Flow Augmentor Fuel Flow 

Ibs/hr. 

WF86 

RALS Augmentor Fuel Flow 

Ibs/hr. 


Symbol 

WR 

W2R 

W22R 

W25R 

X 

XIVEL 

XM 

XM13 

XM93 

XN2 

XN41R 

XN49R 

y 

z 

ZIVEL 

•J?c 

'^1 

9F 

9B 

cr 

30 RF 
60 RF 
80 RF 
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Description 


Engr. Units 


Corrected Air Flow 
Front-Block Fan Corrected Air Flow 
Rear-Block Fan Corrected Air Flow 
Compressor Corrected Air Flow 
State Variable Vector 
Forward Aircraft Velocity 
Flight Mach Number 

Front-Block Fan Discharge Mach Number 
Rear-Block Fan Discharge Mach Number 
Fan Rotor Speed 

Corrected High Pressure Turbine Speed 

m 

Corrected Low Pressure Turbine Speed 

Measurement Error Vector 

Measurement Vector 

Normal Aircraft Velocity 

Compressor Efficiency 

Turbine Efficiency 

ADEN Deflection Angle (Longitudinal) 

RALS Deflection Angle (Longitudinal) 

Standard Deviation 

305» RALS Flow (HTV Transition) 

RALS Flow (HTV Transition) 

80?^ RALS Flow (HTV Transition) 


Ibs/hr. 

Ibs/hr. 

Ibs/hr. 

Ibs/hr. 


t. 

meters/sec. 


RPM 

RPM//^ 


meters/sec. 


Subscripts 


CLEAN 

COMP 

D 

DESIGN 

DIS 

FAN 

IN 


Clean Flow With No Inlet Distortion 
Compressor 
Demand Values 
Nominal Design Value 
Flow With Inlet Distortion 
Fan 

Fan Inlet 
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Syfnbol Dfscriptlon Engr. Units 

OUT Fan Discharge 

Superscripts 

A Most Probable E$t1nMite 

-1 Matrix Inverse 

T Matrix Transverse 
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